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Section  1. 


Executive  Summary 


This  paper  expands  on  the  model  of  sulfur  mustard  (HD)  exposures  and  medical 
countenneasures  (MCMs)  presented  to  the  Defense  Reduction  Agency  (DTRA)  in  201 1, 
(Rodriguez,  Reeves,  &  McClellan,  2011)  to  a  generic  vesicant  model  (GVM).  The  GVM  was 
developed  to  predict  the  effect  of  inhalation,  ocular,  and  percutaneous  exposures  in  healthy 
individuals  and  the  efficacy  of  available  and  proposed  treatments.  The  baseline  models  of 
vesicant  exposure  use  a  series  of  rate  equations  to  describe  the  following: 

•  The  distribution  of  vesicants  at  the  injury  site  and  throughout  the  body. 

•  The  fragmentation  of  DNA  and  suppression  of  glutathione  caused  by  the  exposure  to 
vesicants  in  the  nuclei  of  exposed  cells. 

•  The  series  of  reactions  that  result  in  necrotic  and  apoptotic  epithelial  cell  death  at  three 
different  exposure  sites. 

•  Systemic  oxidative  stress. 

The  MCM  models  modify  the  injury  equations  and  capture  the  therapeutic  effects  of  the 
countenneasure.  In  expanding  the  HD  model  to  the  GVM  we  added  parameters  for  four 
additional  agents,  including  three  nitrogen  mustard  analogues  (HN-1,  HN-2,  and  HN-3),  and 
lewisite.  The  MCMs  included  in  the  framework  are  the  same  as  those  considered  for  the  HD 
model:  bioscavengers,  cell-cycle  inhibitors,  poly(ADP-ribose)  polymerase- 1  (PARP)  inhibitors, 
protease  inhibitors,  and  anti-inflammatories.  We  have  removed  calcium  modulators  from  the 
framework,  as  continued  research  shows  that  the  countenneasures  offer  no  therapeutic  effect. 
Additionally,  we  included  the  chelating  compound  British  anti-lewisite  (BAL),  an  antidote 
against  systemic  lewisite  poisoning. 
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Section  2. 


Introduction  and  Purpose 

Accurate  modeling  of  medical  countenneasure  efficacy  against  chemical,  biological,  and 
radiological  (CBR)  agents  is  essential  to  understanding  the  vulnerability  of  our  warfighters  on 
the  modern  battlefield.  In  helping  calculate  the  benefit  of  countermeasures,  modeling  can  inform 
data-driven  purchasing  decisions  and  logistical  tradeoffs.  In  this  study,  Gryphon  Scientific  and 
Applied  Research  Associates,  Inc.  (ARA)  have  developed  models  to  predict  the  efficacy  of 
medical  countermeasures  against  a  variety  of  agents. 

This  report  (prepared  by  ARA)  is  one  of  ten,  describing  the  medical  countermeasure  models 
constructed  for  this  project.  This  volume  focuses  exclusively  on  the  methods  used  to  construct 
the  Generic  Vesicant  Model  (GVM),  outlines  how  to  use  the  model,  and  provides  examples  of 
the  outputs  generated  by  the  model.  Other  volumes  describe  models  for  B.  anthracis  (volume  1), 
organophosphates  (volume  2),  cesium- 137  (volume  3),  F.  tularensis  (volume  4),  americium-241 
(volume  6),  botulinum  toxin  (volume  8),  plutonium-238/239  (volume  9),  and  plague  (volume 
11).  Each  volume  begins  by  briefly  introducing  the  agent  modeled  and  the  countermeasures 
available  for  use  against  the  agent.  The  overall  schematic  of  each  model  and  the  relevant 
parameters  are  then  discussed,  along  with  a  brief  explanation  of  the  rationale  for  selecting  each 
parameter.  Lastly,  this  report  discusses  the  calculations  and  computational  framework  of  the 
Berkeley  Madonna  model  and  provides  examples  of  modeling  outputs.  We  published  the 
Berkeley  Madonna  source  code  in  the  appendix. 

2 . 1  Summary  of  Deliverables 

Below  is  a  description  of  the  four  deliverables  assigned  for  each  agent  of  interest,  and  a 
description  of  what  is  included  in  each  deliverable.  The  third  numbered  item  ( the  “MCM 
Model”  deliverable)  is  the  deliverable  presented  in  this  interim  report. 

2.  Modeling  Approach 

The  modeling  approach  deliverable  describes  each  of  the  parameters  that  we  anticipate 
including  in  our  MCM  model.  For  each  parameter,  a  description  of  the  approach  for  developing 
and  justifying  the  parameter  is  presented.  The  approach  developed  is  based  on  prior  knowledge 
of  the  agent  and  on  the  general  types  of  data  available  for  each  agent,  but  specific  citations  are 
not  included  as  it  is  a  preliminary  document. 

3.  Modeling  Parameters 

The  modeling  parameters  deliverable  defines  the  value  or  function  for  each  parameter  used  to 
develop  the  MCM  model.  Each  parameter  is  supported  with  a  description  of  the  rationale  for 
choosing  the  parameter,  including  any  scientific  evidence  used  in  parameter  development  or 
assumptions  that  were  made. 

4.  MCM  Model  (This  Report) 

In  addition  to  the  information  already  developed  in  the  “Modeling  Parameters”  deliverable, 
the  “MCM  Model”  includes  a  description  of  the  model,  user  inputs,  the  model  calculations,  and 
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the  model  outputs.  For  the  chemical  agents,  the  report  is  an  accompaniment  to  a  preliminary 
implementation  of  the  MCM  model  built  in  Berkeley  Madonna. 

5.  MCM  Model  in  a  Graphical  User  Interface  (GUI) 

The  MCM  models  will  be  implemented  in  their  final  version  in  a  graphical  user  interface. 
This  implementation  of  the  model  will  include  feedback  and  adjustments  made  after  review  of 
the  previous  deliverables. 

2.2  GVM  Description 

The  basis  for  casualty  and  fatality,  estimation  resulting  from  vesicant  injuries,  will  be  the 
NATO  Planning  Guide  for  Estimation  of  CBRN  Casualties  Ratification  Draft  1  (AMedP-8  (C) 
RD.l)  (NATO,  February  2010).  AMedP-8  (C)  defines  dose  and  time-dependent  sign/symptom 
profiles  for  ocular,  percutaneous,  and  inhalation  exposures  to  HD.  Since  we  are  expanding  this 
model  to  four  additional  vesicants  and  including  the  effects  of  systemic  distribution  and 
oxidative  stress  caused  by  glutathione  suppression,  we  need  to  expand  our  criteria  for  injury  and 
fatalities  beyond  the  sign/symptom  profiles  and  methodologies  defined  in  AMedP-8  (C).  In  this 
document,  we  will  describe  our  proposed  mechanisms  for  incorporating  these  agents  and  injury 
mechanisms,  while  keeping  with  the  thrust  and  goals  of  AMedP-8  (C).  We  have  compared  the 
toxicity  and  injury  mechanisms  of  these  new  agents  with  the  toxicity  and  injury  mechanisms  of 
HD,  to  make  as  much  use  as  possible  of  the  AMedP-8  (C)  injury  profiles.  For  some  agents,  such 
as  the  nitrogen  mustard  analogues,  the  near  equivalence  of  toxicity  and  effect  to  that  of  HD  was  a 
simple  task,  when  compared  to  published  toxicity  data  in  Army  Field  Manual  3.11-9  (FM  3.11- 
9)  (Anny,  2005)  and  the  EPA  (EPA,  Acute  Exposure  Guideline  Levels  for  Selected  Airborne 
Chemicals:  Volume  3,  2003)  (EPA,  2007). 

This  model  represents  a  significant  departure  from  the  HD  model.  We  revisited  the  HD 
framework  and  made  an  almost  complete  change  to  the  skin  injury  model  in  order  to  develop  the 
GVM.  We  did  this  because  we  quickly  realized  that  to  use  the  HD  framework  for  the  GVM,  we 
would  need  to  fit  all  13  rate  parameters,  developed  for  the  HD  model,  for  each  individual  agent. 
However,  most  of  the  rates  that  we  were  modeling  are  completely  independent  of  the  vesicant. 
For  example,  the  rate  of  mitosis  is  a  mechanism  intrinsically  tied  to  the  workings  of  a  cell  and 
will  not  change  when  a  vesicant  is  introduced  into  the  system.  We  revisited  the  literature  we  used 
to  develop  the  HD  model  and  included  some  new  studies.  These  studies  define  equations  that 
represent  the  underlying  biological  mechanisms  in  such  a  way  that  they  are  no  longer  correlated 
with  the  insulting  agent  and,  instead,  are  correlated  with  generic  cellular  injury  and  response. 

The  resulting  model  is  much  more  elegant,  reducing  the  fits  to  three  vesicant-related  rate 
parameters,  only  one  of  which  (the  rate  at  which  the  vesicant  binds  with  and  fragments  the  DNA) 
changes  for  each  insulting  agent.  The  GVM  now  has  a  much  stronger  foundation  in  the  published 
data  and  a  well-defined  process  for  modifying  parameters  as  the  data  becomes  available. 
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Section  3. 


Development  Environment 

The  MCM  model  for  vesicant  exposure  was  developed  and  tested  in  Berkeley  Madonna 
v8.3.18.  Berkeley  Madonna  is  a  differential  equation  solver  ideal  for  modeling  rate  equations, 
allowing  users  to  modify  inputs,  and  perform  batch  runs  to  test  model  sensitivity  and  response. 
Since  the  vesicant  MCM  model  is  a  series  of  rate  equations  designed  to  model  DNA 
fragmentation  and  the  resulting  effects,  Berkeley  Madonna  provides  a  robust  environment  for 
model  development  and  testing. 

Berkeley  Madonna  does  not  export  source  code  into  libraries  that  are  useable  in  C++,  Java, 
or  other  production  languages.  Additionally,  external  programs  cannot  link  into  Berkeley 
Madonna’s  runtime  environment.  We  have  developed  a  C#  differential  equation  solver  that  will 
act  as  a  production  code  framework  for  the  finalized  model.  For  simplicity  of  reproduction  and 
equation  analysis,  we  have  included  the  Berkeley  Madonna  source  code  at  the  end  of  this 
document.  Berkeley  Madonna  is  free  to  download,  but  you  cannot  save  the  source  code  without 
purchasing  a  copy.  However,  you  can  copy  and  paste  the  source  code  into  Berkeley  Madonna  in 
order  to  use  the  model. 

On  January  15th,  2013,  ARA  delivered  the  first  version  of  the  C#  implementation  of  this 
model  to  DTRA  for  review.  On  February  27th,  2013,  a  revised  version  was  delivered  to  DTRA, 
IDA,  and  Gryphon  Scientific.  The  C#  implementation  is  designed  to  run  on  any  Windows  PC 
and  does  not  require  installation. 
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Section  4. 


Vesicant  Injury  Mechanics 

The  model  framework  for  the  GVM  begins  with  normal  cellular  function  in  the  skin  and 
systemic  system  before  vesicant  exposure.  The  model  then  estimates  changes  in  cellular 
function,  in  the  event  of  vesicant  exposure.  The  initial  version  of  this  paper  included  injury  to  the 
eyes  and  the  lungs  in  the  model,  but  ARA  agreed  to  remove  those  systems  from  the  final  model 
after  independent  review  from  the  Institute  of  Defense  Analysis  (IDA),  due  to  a  lack  of 
supportive  data  to  incorporate  these  models.  We  have  preserved  the  discussion  in  Appendix  A. 

4.1  Normal  Cellular  Function 

The  model  is  based  on  a  cell  population  with  four  cellular  states,  summarized  in  Table  4-1. 


Table  4-1.  The  Four  Cellular  States  Included  in  the  Model 


Nn 

Proportion  of  cells  with  uninjured  DNA 

N; 

Proportion  of  cells  with  injured  DNA 

Nk,a 

Proportion  of  killed,  apoptotic  (irreparably  damaged  DNA)  cells 

Nk,n 

Proportion  of  killed,  necrotic  (irreparably  damaged  DNA)  cells 

This  is  a  departure  from  the  HD  model,  which  lumped  apoptotic  and  necrotic  cells  into  a 
generic  “killed  cells”  compartment.  The  differentiation  between  apoptotic  and  necrotic  cells  is 
important  for  two  reasons: 

1.  Necrotic  cells  result  from  the  severe  depletion  of  NAD+  and  ATP,  which  in  turn  signals 
protease  and  leads  to  vesication  and  severe  blistering  (Kehe,  Balszuweit,  Steinritz,  & 
Thiermann,  2009)  (Kehe  &  Szinicz,  2005).  Since  severity  of  vesication  is  the  metric  for 
injury  severity  in  our  model,  this  is  an  important  distinction.  Whereas  we  won’t  see  large 
differences  in  untreated  individuals  following  vesicant  exposure,  we  will  see  large 
differences  in  individuals  treated  with  PARP  inhibitors. 

2.  PARP  inhibitors,  which  have  been  shown  to  have  no  effect  on  the  gross  viability  of  cells  or 
level  of  vesication  (Atkins,  Lodhi,  Hurley,  &  Hinshaw,  2000),  but  have  a  great  effect  on  the 
distribution  of  necrotic  and  apoptotic  cells  (Kehe  K.  ,  2008).  In  the  HD  model 
implementation,  we  saw  that  our  model  predicted  no  therapeutic  effect  when  PARP 
inhibitors  were  administered.  However,  for  skin  injuries,  the  therapeutic  effect  can  actually 
be  significant  if  vesication  is  minimized. 

For  this  implementation  of  the  model,  we  only  consider  percutaneous  injuries  and  systemic 
injuries.  For  systemic  injuries,  we  will  only  consider  the  systemic  concentration  of  vesicants  in 
this  model  and  correlate  the  concentration  of  vesicant  to  acute  lethality  using  probit  analysis 
from  the  toxicity  data  presented  in  FM  3.11-9  (Army,  2005) 
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We  will  routinely  use  the  tenn  “relative  concentration”  throughout  this  report.  Essentially, 
the  model  does  not  attempt  to  estimate  the  actual  molar  concentration  of  any  chemical,  molecule, 
or  cell  included  in  the  framework.  Instead,  it  attempts  to  measure  the  increase  (or  decrease)  of 
concentrations  for  each  chemical,  molecule,  or  cell  included  in  the  model.  In  order  to  achieve 
this,  we  set  initial  conditions,  so  that  any  chemical,  model,  or  cell  that  increases  in  concentration 
as  a  result  of  vesicant  exposure  has  an  initial  relative  concentration  of  0.0.  Any  chemical,  model, 
or  cell  that  decreases  as  a  result  of  vesicant  exposure  has  an  initial  relative  concentration  of  1.0. 

The  model  assumes  that,  prior  to  vesicant  exposure,  the  relative  concentrations  of  normal, 
healthy  cells  in  the  skin  are  1.0  and  the  relative  concentration  of  injured  and  killed  cells  is  0.0. 
These  values  signify  that  100%  of  the  cells  in  the  area  of  interest  are  nonnal,  0%  are  injured,  and 
0%  are  killed.  The  model  does  not  attempt  to  estimate  actual  cell  concentrations,  rather  the 
fraction  of  cells  above  or  below  the  normal  percentage.  In  the  event  of  injury,  the  relative 
concentration  of  nonnal  cells  will  decrease  while  fractured  and  killed  cells  will  increase.  During 
the  recovery  phase,  the  relative  concentration  of  normal,  healthy  cells  will  approach  1,  while  the 
relative  concentration  of  injured  cells  will  approach  0.  The  relative  concentration  of  killed  cells 
is  cumulative  and  will  never  return  to  1 . 

In  addition  to  the  three  cell  states,  the  model  also  predicts  relative  concentrations  of  PARP, 
NAD,  protease,  inflammatory  cells,  glutathione  (GSH),  ATP,  and  calcium  ions.  All  of  these 
chemicals  are  assumed  to  be  present  in  each  cell  before  vesicant  exposure,  making  their  relative 
concentrations  equal  to  either  1.0  or  0.0,  depending  on  whether  their  concentration  increases  or 
decreases.  Following  exposure,  some  of  these  elements  will  decrease  to  zero  (NAD,  for  instance, 
which  depletes  during  DNA  repair)  and  some  of  these  elements  will  increase  to  1.0  or  higher 
(protease,  for  instance,  which  is  used  to  break  down  killed  cells). 

4.2  Injury  Mechanisms 

The  differences  in  vesicant-related  skin,  eye,  lung,  and  other  tissue  injuries  were  discussed  in 
the  HD  model.  We  will  not  discuss  the  ocular  and  inhalation  injury  mechanisms  in  this  paper,  as 
they  are  not  included  in  the  final  GVM.  We  will  discuss  enhancements  to  the  skin  and  systemic 
systems;  however,  that  is  meant  to  augment  the  discussion  in  the  original  HD  paper.  The  types  of 
injured  and  killed  cells  considered  in  the  model  are  listed  in  Table  4-2. 


Table  4-2.  Types  of  Injured  and  Killed  Cells 


Variable 

Cell  Type 

Nr 

Injured  skin  cell 

Nka 

Apoptotic  skin  cell 

Nk.n 

Necrotic  skin  cell 

For  systemic  injuries,  we  will  only  consider  the  concentration  of  vesicant  in  the  circulatory 
system  in  this  model  and  correlate  the  concentration  to  acute  lethality,  described  later. 
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4.2.1.  Skin 


The  skin  and  mucous  membranes,  which  cover  most  of  the  body  and  account  for  over  10%  of 
its  mass,  provide  a  protective  layer  against  environmental  toxins,  injurious  substances,  and 
pathogens.  The  outer  layer  of  the  skin  (stratum  comeum)  functions  as  a  barrier  and  prevents  most 
materials  from  absorption  through  the  lower  epidennal  layers  into  the  dermis  and  deeper  tissues, 
including  the  systemic  circulation  (Centers  for  Disease  Control  and  Prevention  (CDC),  2010). 
Figure  4-1  shows  the  different  layers  of  the  skin.  In  the  GVM,  we  now  account  for  the 
distribution  of  vesicant  across  these  different  layers,  focusing  on  the  stratum  corneum,  the  basal 
epidermis,  the  dermis,  and  the  vascular  system. 
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Figure  4-1.  Cross-Section  of  Skin  (SEER  Training:  Layers  of  the  Skin) 

As  mentioned  in  the  original  FID  paper,  vesication  and  severe  blistering  are  primarily 
achieved  when  vesicants  enter  the  basal  epidermis  and  dermis,  alkalizing  the  healthy  DNA  and 
causing  necrotic  cell  death.  The  toxicokinetic  model  that  we  implemented  for  the  GVM  and 
describe  in  section  4.3  estimates  the  proportion  of  HD  applied  to  the  skin  that  ends  up  in  the 
basal  epidermis  and  the  dermis,  leading  to  vesication.  Additionally,  we  included  a  mechanism 
that  accounts  for  GSH  suppression,  which  has  been  shown  to  increase  oxidative  stress  and  acute 
toxicity.  We  made  another  enhancement  to  the  injury  mechanism  in  the  skin  to  differentiate 
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between  apoptotic  and  necrotic  cell  death,  for  reasons  stated  in  the  preceding  section.  Once 
again,  we  make  this  distinction  in  the  basal  epidermis  and  dermis  layers  of  the  skin  and  correlate 
necrotic  cell  death  to  vesication. 

The  mechanism  that  differentiates  apoptotic  and  necrotic  cell  death  is  still  being  studied,  but 
sufficient  new  information  was  published  to  allow  us  to  include  the  differentiation  into  the 
model.  We  found  several  earlier  studies  that  looked  at  the  effect  of  N-acetylcysteine  (NAC, 
which  is  included  in  our  model  as  a  bioscavenger)  against  both  nitrogen  mustard  (Weltin,  et  ah, 
1996)  and  sulfur  mustard  (Greenberg,  Kamath,  Petrali,  Hamilton,  Garfield,  &  Garlick,  2006) 
(Atkins,  Lodhi,  Hurley,  &  Hinshaw,  2000)  injuries  with  data  that  originally  led  us  to  believe  that 
apoptosis  was  the  primary  injury  mechanism  for  vesication  and  blistering. 

However,  more  recent  papers  have  shown  that  necrosis  plays  a  very  important  role  in 
vesication  and  is  the  cause  of  severe  injuries  in  people  exposed  to  vesicants.  Nicotera  &  Melino 
(Nicotera  &  Melino,  2004),  who  were  not  working  with  vesicant  exposures,  presented  evidence 
that  the  apoptotic-necrotic  switch  was  a  complicated  process,  originally  believed  to  be  tied 
exclusively  to  capsases  (the  HD  model  we  developed  even  used  protease  concentrations  as  a 
precursor  to  cell  death  and  vesication),  but  now  shows  significant  evidence  that  it  is,  in  fact,  tied 
to  ATP  depletion.  The  paper  presents  a  complicated  model  of  ATP  regulation,  but  for  our 
purposes,  we  can  simply  focus  on  relative  ATP  concentration  in  order  to  shift  the  weighting 
between  apoptotic  and  necrotic  cell  death.  Figure  4-2,  which  was  taken  from  Nicotera  & 
Melino’s  paper,  shows  the  simplified  model  of  apoptosis/necrosis  dependence  on  ATP  levels. 
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ATP-dependentPS  franslocase 
/ipoptosome  activation 
A'I'm  stability 
Calpain  activation 
Other  targets 


Necrosis 


•  Random  morphology 

•  Passive  execution 
•Caspases  not  necessarily  involved 

•  Frequent  in  pathological  settings 

>.  -wm  t: 


Apoptosis 

•Typical  morphology 

•  Active  execution 
•Caspase  execution 
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Figure  4-2.  Apoptosis-Necrosis  Switch  by  ATP  (Nicotera  &  Melino,  2004) 
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Tying  this  into  the  GVM,  our  model  accounts  for  NAD+  depletion  which,  in  turn,  leads  to  a 
depletion  of  ATP.  We  did  not  include  ATP  in  the  HD  model,  but  we  now  include  it  in  the  GVM. 
The  vesicant  research  community  has  recently  published  papers  that  also  strengthen  the 
argument  for  ATP  levels  being  a  regulator  of  necrosis  and  apoptosis  and  also  detennines  how  the 
differentiation  relates  to  the  type  of  vesicant  injury.  Figure  4-3  is  from  a  paper  published  by 
Kehe,  et.  al  (Kehe,  Balszuweit,  Steinritz,  &  Thiennann,  2009)  of  a  study  that  aimed  to  show  how 
NAD+/ATP  depletion  and  necrosis  is  concentration-dependent,  meaning  large  doses  of  HD  lead 
to  a  severe  depletion  of  NAD+  and  ATP,  which  in  turn  lead  to  an  increase  in  necrotic  cells. 


Figure  4-3.  Apoptosis  and  Necrosis  Determination  Following  HD  Exposure  (Kehe, 
Balszuweit,  Steinritz,  &  Thiermann,  2009) 

Another  recent  study  (Debiaka,  Keheb,  &  Burkle,  2009)  also  correlated  NAD+  depletion  to 
necrosis  and  blister  fonnation  in  HD  exposures.  Figure  4-4,  which  was  taken  from  the  paper, 
once  again  shows  how  PARP  activity  leads  to  NAD+  depletion,  necrosis,  and  eventually 
inflammation  and  tissue  damage.  In  the  Debiaka,  et.  al.  model,  there  is  an  additional  pathway  to 
necrosis  presented  that  comes  from  imperfect  DNA  repair.  This  pathway  could  be  considered  for 
future  investigation  and  research  in  future  development  of  the  GVM. 
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Figure  4-4.  Detalied  Diagram  of  Apoptosis/Necrosis  Related  to  HD  Injury  (Debiaka,  Keheb, 

&  Burkle,  2009) 


Also,  in  a  recent  study  of  human  exposures  (Pita  &  Vidal -Asensi,  2010),  the  authors  state, 

“In  severe  lesions,  particularly  those  due  to  exposure  to  the  liquid  agent,  necrosis  occurs  and  they 
become  more  susceptible  to  secondary  infections  that  can  complicate  the  course  of  the  lesion.” 
This  statement  strengthens  the  argument  that  large  doses  result  in  blister  fonnation  and, 
eventually,  necrosis.  In  the  same  paper,  the  author  also  states  the  outer-most  layer  of  some 
blisters  in  human  injuries  is  “completely  necrotic.” 

4.2.2.  Systemic 

The  toxicokinetic  model  also  allows  us  to  measure  the  proportion  of  vesicant  that  enters  the 
circulatory  system.  We  are  not  yet  including  a  model  of  systemic  injury,  but  we  are  modeling 
GSH  suppression  in  following  systemic  distribution  that  occurs  primarily  in  the  bone  marrow, 
gastrointestinal  tract,  and  the  CNS  (Pitaa  &  Vidal-Asensib,  2010).  The  bone  marrow  suppression 
leads  to  the  second  cause  of  death  outside  of  acute  respiratory  failure.  Due  to  time  constraints 
and  a  current  lack  of  data,  we  will  not  model  bone  marrow  suppression  and  the  subsequent 
bacterial  pneumonia,  but  will  instead  tie  all  of  systemic  injuries  to  oxidative  stress  caused  by  the 
systemic  dose.  The  model  being  used  here  is  likely  incomplete,  but  the  mechanism  of  correlating 
a  percutaneous  dose  to  systemic  dose  to  lethality  is  well  founded.  Our  first  step  will  be  to  use  the 
toxicokinetic  model  to  estimate  the  quantity  of  applied  dose  that  ends  up  in  the  circulatory 
system.  The  second  step  would  be  to  predict  the  level  of  glutathione  suppression  resulting  from 
this  systemic  dose.  The  third  step  would  be  to  use  the  toxicity  values  published  by  ECBC,  for  all 
of  the  vesicants  of  interest,  to  correlate  systemic  dose  and  glutathione  suppression  levels  to  the 
LD5o  and  probit  for  percutaneous  dose.  Although  we  are  missing  the  finer  resolution  that  a 
complete  systemic  model  would  require,  we  believe  that  we  have  a  model  that  ties  systemic  dose 
to  lethality.  The  model  will  allow  us  to  account  for  chelating  therapy  using  British  anti-lewisite 
(BAL),  which  removes  lewisite  from  the  circulatory  system  and  will  also  allow  us  to  better 
include  the  beneficial  effects  of  bioscavengers  like  NAC,  which  work  to  increase  the  level  of  free 
GSH  in  the  body. 
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Section  5. 


Vesicant  Exposure  Model 

Unlike  with  the  HD  model,  we  are  only  considering  percutaneous  exposures  in  the  GVM.  In 
this  section,  we  will  revisit  the  percutaneous  exposure  models  as  we  made  a  change  to  the 
original  HD  model  to  account  for  a  correction  in  the  MOPP  calculation  and  we  also  included  a 
toxicokinetic  model  for  percutaneous  distribution  in  the  skin  and  the  blood  stream.  We  have 
preserved  the  Inhalation  and  Ocular  exposures  in  Appendix  A. 

5.1  Percutaneous 

For  the  GVM,  we  have  expanded  the  HD  percutaneous  model  to  include  a  toxicokinetic 
model  that  predicts  vesicant  concentration  at  different  stratum  of  the  skin  and  vesicant  that  enters 
the  circulatory  system  and  is  distributed  systemically.  This  is  an  important  feature  for  vesicants, 
as  the  delayed  distribution  time  of  the  chemical  can  help  the  model  better  predict  delayed  onset 
of  effects.  The  percutaneous  model  will  now  make  use  of  a  compartmental  model  that  was 
developed  by  Cutaneous  Phannacology  and  Toxicology  Center  at  North  Carolina  State 
University  (Riviere  &  al.,  1994),  illustrated  in  Figure  5-1. 
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Figure  5-1.  Multi-Compartmental  Toxicokinetic  Model  for  HD  Penetration  in  the  IPPSF 

(Riviera  &  al.,  1994) 

The  model  was  developed  using  the  isolated  perfused  porcine  skin  flap  (IPPSF)  as  an  analog 
for  human  skin.  IPPSF  was  used  to  study  the  percutaneous  absorption  of  many  different  drugs 
and  chemicals  (Riviere  &  Monteiro-Riviere,  1991)  (Williams,  Carver,  &  Riviere,  1990).  One 
experiment  that  documented  the  absorption  of  HD  using  the  IPPSF  was  done  by  the  Cutaneous 
Pharmacology  and  Toxicology  Center  at  North  Carolina  State  University  in  1997  (Monteiro- 
Riviere  &  Inman,  1997).  In  that  study,  the  morphological  similarities  between  human  injuries 
and  HD-induced  injury  on  the  IPPSF  are  very  similar.  IPPSF  is  the  only  human-skin  analogue 
found,  to  date,  that  actually  produced  gross  blisters  when  exposed  to  HD,  much  like  human 
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injuries.  The  model  shown  in  Figure  5-1  has  a  series  of  rates  that  have  been  derived  from 
experimental  data  and  are  given  in  a  Riviere,  et  al.  study  (Riviera  &  ah,  1994).  The  rates  are 
reproduced  in  Table  5-1. 


Table  5-1.  Rate  Constants  for  IPPSF  HD  Deposition  Model  (Riviera  &  al.,  1994) 


Constant 

1<48 

]<49 

k-94 

kq2 

k29 

k23 

k?2 

a 

b 

ki3 

k35 

k-53 

K-2,22 

K-22,2 

Value 

(hr) 

0.34 

0.25 

0.23 

0.26 

0.22 

0.24 

0.12 

0.33 

0.0001 

0.12 

0.0058 

0.006 

0.12 

0.007 

The  rate  constants  are  used  in  a  series  of  equations,  shown  below,  that  predict  the  dose  of 
vesicant  at  each  of  the  compartments.  These  equations  were  not  defined  in  the  Riviera  study,  but 
they  take  the  form  of  standard  rate  equations  for  toxicodynamic  processes  and  we  verified  the 
model  output  (Figure  5-2  -  Figure  5-6)  against  the  published  results.  The  equations  are  shown 
below. 


The  concentration  of  evaporated  mustard  ( Evap )  off  of  the  surface  of  the  skin  (Surf): 


d 

dt 


Evap  =  k48  *  Surf 


(5-1) 


The  concentration  of  mustard  on  the  surface  of  the  skin  and  the  relationship  between  the 
stratum  corneum  concentration  (SCorn): 


d 

—Surf  =  k94  *  SCorn  —  k48  *  Surf—k49  *  Surf 
dt  ‘  L 


(5-2) 


The  concentration  of  mustard  in  the  stratum  comeum  and  the  relationship  between  the 
surface  and  the  basal  epidermis  concentration  ( BEpi ): 

d  (5-3) 

—SCorn  =  k49  *  Surf  +  k29  *  BEpi  —  k92  *  SCorn  —  k94  *  SCorn 
dt 

The  concentration  of  mustard  in  the  basal  epidermis  and  the  relationship  between  the  stratum 
comeum,  dennis  (Derm)  and  metabolites  (Met)  concentrations: 

d  (5-4) 

—  BEpi  =  k92  *  SCorn+k32  *  Derm+k22  2  *  Met  —  k23  *  BEpi—k29  *  BEpi—k2  22  *  BEpi 
dt 


The  concentration  of  mustard  in  the  dermis  and  the  relationship  between  the  basal  epidermis, 
vascular  system  (Vas),  and  fat  (Fat)  concentrations: 

d  h  ,  <5-5> 

—  Derm  =  k23  *  BEpi  +  k13  *  Vas  +  k53  *  Fat  —  a  *  e  *  Derm  —  k32  *  Derm  —  k35 
dt 

*  Derm 
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The  concentration  of  mustard  in  the  fat  and  the  relationship  between  the  dennis 
concentration: 


d  (5-6) 

~rFat  =  k3S  *  Derm  —  k53  *  Fat 


The  concentration  of  mustard  in  the  metabolites  and  the  relationship  between  the  basal 
dennis  concentration: 


dt 


Met  —  /c 2  22  *  EEpi  k32  2  *  Met 


(5-7) 


The  concentration  of  mustard  in  the  vascular  system  and  the  relationship  between  the  dermis 
concentration: 


ci 

—  Vas  =  a*  e  *  Derm—k13  *  Vas 


(5-8) 


dt 


This  just  leaves  an  equation  for  effluent  mass.  We  turned  to  a  more  recent  paper  (Riviere  & 
al,  2001)  for  this  calculation  and  made  use  of  the  model: 


T(t)=4(e  t/zb b)— e  f/Td) 


(5-9) 


Where: 

1.  Y(t)  is  an  estimate  of  the  effluent  mass.  This  will  be  known  as  the  systemic  dose,  Ds. 

2.  ^4  is  the  absorbed  dose  or  the  mass  in  the  dennis, 

3.  The  first  exponential  term  represents  the  rate  of  tenninal  release  of  the  absorbed  chemical, 
and 

4.  The  second  exponential  term  represents  the  rate  of  uptake  of  the  chemical  from  the  skin  to 
the  IPPSF 

The  values  of  Xb  and  id  were  obtained  from  Riviera,  2001,  and  are  equal  to  39  minutes  and 
0.7  minutes,  respectively.  We  ran  the  model  with  an  input  surface  dose  of  3000  pg  of  HD,  as  was 
done  in  Riviera,  1994,  in  order  to  gauge  the  model  output  compared  to  the  calculated  output.  The 
result  of  our  output  is  in  Figure  5-2  through  Figure  5-6.  When  compared  to  the  output  from  the 
referent,  we  see  a  near  perfect  reproduction  of  the  model  output  (which,  in  turn,  reproduces  the 
experimental  data  as  shown  in  the  referent)  for  all  compartments,  except  for  the  effluent  mass, 
which  is  higher  in  our  model  (final  mass  at  8  hours  is  equal  to  180  pg  instead  of  1 10  pg).  We 
attribute  this  to  using  an  updated  effluence  model  for  our  implementation,  as  mentioned  above. 
The  next  section  documents  how  we  use  these  calculated  exposures  in  the  mustard  injury  model. 
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Calculated  Stratum  Corneum  and  Total  Skin  Mass 


Time  (Hours) 

Figure  5-3.  Calculated  Stratum  Corneum  and  Total  Skin  Mass  for  3000jig  of  HD 
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Calculated  Dermis  and  Effluent  Mass 


Time  (Hours) 


Figure  5-4.  Calculated  Dermis  and  Effluent  Mass  for  3000jag  of  HD 
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Calculated  Basal  Epidermis  and  Metabolite  Mass 


Time  (Hours) 

Figure  5-6.  Calculated  Basal  Epidermis  and  Metabolite  Mass  for  3000pg  of  HD 

5.2  Using  the  Calculated  Dose  in  the  GVM 

For  the  percutaneous  and  systemic  dose,  using  the  calculated  dose  modified  the  HD  model’s 
implementation  significantly.  These  changes  are  discussed  in  this  section. 

One  note  on  dose:  we  assess  the  dose  and  the  effect  at  each  injury  site  separately.  There  are 
two  doses  we  calculate,  summarized  in  Table  5-2. 


Table  5-2.  The  Two  Doses  Included  in  the  GYM 


Dose 

Ds 

DP 

Value 

Systemic  Dose 

Percutaneous  Dose 

The  systemic  and  percutaneous  doses  have  two  sink  terms,  whole  DNA  in  nonnal  cells  and 
free  GSH. 

For  percutaneous  dose,  the  first  step  is  to  calculate  the  combined  dose  deposited  on  the  skin. 
This  is  given  by  the  equation  below: 

d  Hmask  rTMOPP  (5-10) 

—  Deposited.  =  C  dt  *  CF  +  I  C  dt  *  CF  *  Sexposed 

^  Tmask 
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This  equation  estimates  the  whole-body  percutaneous  vapor  exposure,  up  until  a  mask  is 
donned  (Tmask).  Plus,  estimates  of  the  percutaneous  vapor  exposure  to  the  proportion  of  the  body 
left  unprotected  (Sexposed=  0.93),  up  until  the  time  at  which  the  rest  of  MOPP  gear  is  donned 
(Tmopp)-  The  other  tenn  used  above  is  the  conversion  factor  from  percutaneous  vapor  to  an 
equivalent  percutaneous  liquid  dose,  CF.  In  the  HD  paper  we  explained  that  the  value  is  equal  to 
the  ratio  of  the  ED50  for  severe  percutaneous  liquid  exposure  to  the  ECt50  for  severe 
percutaneous  vapor  exposure.  This  is  an  agent-specific  value  that  is  summarized  in  Table  5-3. 
Values  for  the  ED5q  and  ECt5o  were  obtained  from  FM  3-1 1.9  (Army,  2005). 


Table  5-3.  Conversion  Factors  for  Vesicants  Included  in  the  GVM 


Agent 

HD 

HN-1,  2,  and  3 

Lewisite 

ED50,  Severe  Percutaneous,  Liquid  (mg/man) 

600 

600 

600 

ECt50,  Severe  Percutaneous,  Vapor  (mg-min/m3) 

500 

500 

500 

CF  (m/min) 

1.2 

1.2 

1.2 

Equation  (  5-10  )  is  initialized  to  be  equal  to  the  liquid  skin  dose  in  milligrams.  The  result  of 
this  sub-model  can  be  seen  in  Figure  5-7,  where  an  individual  was  exposed  to  10  mg  of  liquid 

■J 

vesicant  and  a  constant  concentration  of  1  mg/m  of  vapor  vesicant,  until  a  mask  was  donned  at  2 
minutes  and  full  MOPP  gear  was  donned  at  5  minutes.  As  you  can  see,  the  mask  makes  only  an 
imperceptible  change  in  the  accumulating  dose  for  percutaneous  dosing. 


Percutaneous  Dose  Deposited  on  Skin 
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Figure  5-7.  Equivalent  Percutaneous  Dose  Deposited  on  Skin  for  a  Constant  Vapor 
Exposure  and  Instantaneous  Liquid  Exposure 
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The  amount  of  vesicant  deposited  on  the  skin  (. Deposited)  is  then  used  in  the  toxicokinetic 
algorithm  by  using  it  as  the  initial  value  equation 

( 5-2  ),  Surf. 

Figure  5-8  shows  the  dose  at  the  stratum  comeum,  basal  epidermis,  dermis,  and  vascular 
compartments,  as  well  as  the  total  skin  dose,  as  a  function  of  time  for  the  inputs  used  in  Figure 
5-7,  calculated  in  Berkeley  Madonna. 


Dose  at  Different  Skin  Layers 


Time  (Minutes) 

Figure  5-8.  Dose  Over  Time  at  Four  Different  Skin  Layers 

The  benefit  of  the  toxicokinetic  model  is  its  accounting  for  the  distribution  of  vesicant  across 
the  skin,  allowing  us  to  better  predict  the  extent  of  an  injury  against  a  given  dose.  Vesicants, 
being  lipophilic,  can  diffuse  readily  through  the  intracellular  spaces  of  the  stratum  comeum, 
which  are  filled  with  lipids  such  as  sebum,  oils,  waxes,  as  well  as  sweat  (largely  water).  From 
there  vesicants  can  make  their  way  through  the  intermediate  layers  into  the  basal  epidermis, 
where  the  basal  cells  (which  produce  new  skin  cells,  among  other  functions)  are  located.  Injuring 
the  basal  epidermis  and  the  dermis  is  what  ultimately  results  in  vesication,  which  is  the  metric 
for  skin  injury  in  our  model.  For  this  reason,  the  percutaneous  dose,  Dp,  is  the  sum  of  the  Basal 
Epidermis  and  Dermis  compartment  or: 

Dp(t)  =  BEpi(t )  +  Derm(f)  (  $~H  ) 
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It  is  important  to  note  that,  unlike  the  ocular  and  inhaled  doses,  we  will  not  incorporate  a  sink 
tenn  for  Dp.  Toxicodynamic  models  already  account  for  sink  terms  when  the  mass  of  an  agent  is 
calculated  at  different  compartments.  For  the  compartment  model  used  here,  radiolabeled  HD 
was  administered  and  monitored  at  the  various  skin  layers.  The  dose  was  fragmenting  DNA  or 
suppressing  GSH,  as  it  was  there.  The  dose  will  be  incorporated  into  the  model  by  fragmented 
DNA  and  suppressing  GSH,  while  it  is  present  in  the  basal  epidennis  or  dennis,  but  it  will  not 
persist  to  injure  new,  healthy  cells  until  all  of  the  vesicant  is  used,  as  was  the  case  in  the  HD 
model. 

The  final  dose  is  the  systemic  dose,  Ds.  The  systemic  dose  is  the  effluent  dose  calculated  by 
the  toxicokinetic  model  or: 


Ds(i)  =  Y(t)-kG*G*Ds(i)  (5-12> 

For  systemic  dose,  we  are  primarily  concerned  with  GSH  suppression  and  acute  lethality  via 
the  percutaneous  route;  we  ignored  these  two  factors  in  the  HD  model.  We  will  include  GSH  as  a 
sink  term  for  systemic  dose  and  use  GSH  suppression  as  the  input  to  the  acute  lethality 
calculation  that  will  be  explained  later  in  the  document. 
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Section  6. 


Vesicant  Injury  Model  Framework 

Figure  6-1  shows  the  injury  model  framework  that  was  developed  for  FID.  This  model  was 
used  for  ocular,  inhalation,  and  percutaneous  exposures.  For  the  GVM,  we  are  modifying  this 
model  slightly  for  ocular  and  inhalation  exposures  (Figure  6-1)  and  we  are  expanding  the  injury 
model  framework  for  percutaneous  exposures  to  include  GSH  suppression  for  both  percutaneous 
exposures  and  systemic  injuries  (Figure  6-2). 


Figure  6-1.  Model  Framework  Developed  for  HD 


For  the  percutaneous  model  (Figure  6-2),  we  made  three  significant  changes  to  the 
framework: 

1 .  Differentiating  between  necrotic  and  apoptotic  cellular  death.  This  will  allow  us  to  better 
estimate  the  effect  of  PARP  inhibitors,  which  have  been  shown  to  have  negligible  effects  on 
the  percentage  of  viable  cells  following  mustard  exposure,  but  have  been  shown  (Kehe  K., 
Inhibition  of  polu(ADP-ribose)  polymerase  (PARP)  influences  the  mode  of  sulfur  mustard 
(SM)-induced  cell  death  in  HaCaT  cells,  2008)  to  greatly  change  the  ratio  of  necrotic  and 
apoptotic  deaths,  thus  decreasing  the  severity  of  vesication. 
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2.  The  addition  of  ATP  and  its  depletion  being  the  primary  signaling  pathway  for  protease. 

3.  Systemic  toxicity  that  results  from  the  vesicant  entering  the  bloodstream.  This  will  allow  us 
to  account  for  acute  lethal  percutaneous  doses  of  vesicants  as  suggested  by  FM3-1 1.9. 

Relating  the  compartments  in  the  percutaneous  model  to  the  discussion  in  the  previous 
section  on  exposure,  the  vesicant  liquid  and  vapor  combine  to  form  the  initial  surface  dosage, 
which  either  evaporates  or  absorbs  into  the  skin.  Within  the  compartment  labeled  “Combined 
Skin  Exposure”  is  the  toxicokinetic  model.  The  basal  epidermis  and  dennis  doses  combine  to 
attack  the  DNA,  while  the  effluent  dose  enters  systemic  distribution. 


Figure  6-2.  GVM  Model  Framework  for  Percutaneous  Exposures 


An  additional  benefit  to  the  revised  framework  is  an  easier  discussion  on  the  individual 
components,  as  we  can  now  look  at  each  box  as  subcomponents  of  the  model  and  provide 
differential  equations  for  each  one.  The  remainder  of  this  section  will  present  and  discuss  the 
system  diagrams  and  equations  for  each  subcomponent.  We  will  also  present  the  data  that  was 
used  to  derive  the  parameters  for  the  different  equations. 

6.1  Injury  Model  Overview 

We  once  again  have  a  significant  departure  from  the  HD  model.  In  the  HD  model,  PARP, 
NAD+,  protease,  Ca++,  and  inflammatory  cells  were  modeled  the  same  way  as  normal,  injured, 
and  killed  cells  -  meaning  they  had  source  and  sink  terms  and  various  differential  equations  to 
describe  their  relative  concentrations.  Revisiting  the  literature  and  the  framework  allowed  us  to 
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better  understand  the  complex  relationships  between  these  molecules  and  lead  to  a  better 
understanding  of  how  they  work.  We  believed  we  were  modeling  them  incorrectly  after  further 
review. 

For  example,  we  modeled  PARP  as  a  differential  equation  that  was  dependent  on  two  tenns: 

1 .  PARP  could  be  added  to  the  system  at  a  rate  kPARPrel,  dependent  on  the  concentration  of 
injured  cells. 

2.  PARP  could  be  removed  from  the  system  as  it  repairs  damaged  DNA  at  a  rate  krep, 
dependent  on  the  concentration  of  injured  cells  and  NAD+ 

There  is  a  fallacy  in  this  equation,  as  PARP  is  never  truly  “removed”  from  the  system, 
instead  it  is  momentarily  used  to  repair  damaged  DNA.  One  solution  to  this  problem  would  be  to 
introduce  a  third  term  that  accounts  for  PARP’s  reintroduction  into  the  system;  however,  there  is 
still  no  way  to  return  PARP  back  to  its  baseline  value  of  0.0,  unless  we  introduce  a  fourth  term 
that  moves  PARP  back  to  its  baseline  concentration.  Using  Berkeley  Madonna  to  fit  two 
unknown  rates  to  one  data  set  was  a  difficult  task  to  justify,  but  fitting  four  rates  to  one  data  set 
would  not  result  in  a  good  model. 

Our  solution  to  this  problem  was  to  once  again  look  at  the  framework  and  the  available  data 
and  state  several  assumptions  to  better  construct  our  equations: 

1 .  PARP  concentration  is  only  dependent  on  the  concentration  of  injured  cells 

2.  NAD+  concentration  is  only  dependent  on  the  concentration  of  PARP 

3.  ATP  concentration  is  only  dependent  on  the  concentration  of  NAD+ 

4.  Protease  is  only  dependent  on  the  concentration  of  ATP 

5.  Inflammatory  cells  and  calcium  ions  are  only  dependent  on  the  concentration  of  killed  cells 

This  may  seem  like  a  simplification  of  the  model,  but  it  is  consistent  with  the  injury 
mechanism,  in  turn,  consistent  with  the  data,  and  reduces  the  fits  down  to  nine  equations  with 
three  unknowns.  Additionally,  it  keeps  with  the  philosophy  of  the  model:  that  by  incorporating 
the  biochemical  processes  that  lead  to  injury,  we  can  better  study  how  to  disrupt  those  processes 
with  countermeasures.  The  new  framework  implementation  is  built  from  the  ground  up,  starting 
with  the  biochemical  processes  (which  are  independent  of  agent  and  exposure)  and  overlaying 
the  agent  and  exposure-dependent  equations  on  top  of  them. 

We  now  present  the  equations  and  parameters  for  the  underlying  biological  processes  and  the 
cellular  response  to  these  processes. 

6.2  Underlying  Biological  Processes 

This  section  details  the  equations  and  data  that  were  used  to  provide  fits  for  the  underlying 
biological  processes:  PARP,  protease,  inflammatory  cells,  calcium  ion  signaling,  and  NAD+  and 
ATP  depletion.  We  assume  that  these  processes  are  agent-independent. 

6.2.1.  PARP,  NAD+,  and  ATP  Concentrations 

In  order  to  define  equations  for  PARP,  NAD+,  and  ATP  we  needed  data  correlating  injured 
cells  to  each  of  these  parameters.  Whereas  there  are  many  studies  that  examine  these  molecules 
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independently,  we  used  one  study  that  looked  at  cell  injuries  resulting  from  HD  exposure  and  the 
effect  these  injuries  have  on  PARP,  NAD+,  and  ATP  levels.  This  study,  performed  by  Daniel 
Hinshaw,  et  al.  at  the  Department  of  Veterans  Affairs  Medical  Center  (Hinshaw,  Lodhi,  Hurley, 
Atkins,  &  Dabrowska,  1999),  looked  at  in  vitro  exposures  to  human  endothelial  cells  and 
keratinocytes,  to  better  understand  the  role  that  PARP  plays  in  the  injury  process.  In  doing  so, 
the  researchers  managed  to  look  at  the  extent  of  injury  to  the  cells  and  the  PARP,  NAD+,  and 
ATP  levels  as  a  function  of  time  and  dose.  Whereas  this  study  isn’t  very  recent  (  published, 
1999),  it  best  represents  the  type  of  study  that  can  be  used  to  provide  data  for  the  GVM.  We 
decided  to  use  the  study  for  this  parameter  fit  as  the  data  is  complete  and  allows  us  to  define  the 
process  and  the  framework  for  populating  our  model. 

As  mentioned,  the  model  studied  both  endothelial  cells  and  keratinocytes.  We  focused  on  the 
keratinocytes  since  they  make  up  the  majority  of  the  cells  found  in  the  basal  epidennis  and 
dennis  layers  that  we  are  interested  in.  The  study  looked  at  extent  of  injury  and  PARP,  NAD+, 
and  ATP  levels  against  0,  250,  500,  and  1000  pM  of  HD  at  2,  3,  5,  and  6  hours  post-exposure. 
Unfortunately,  there  weren’t  perfect  intersections  across  all  time-points  -  for  instance,  cell  injury 
was  not  measured  at  2-hours  and  PARP,  NAD+,  and  ATP  were  not  measured  at  6-hours,  but  we 
were  still  left  with  enough  data  to  develop  a  process  and  initial  fits. 

Additionally,  the  study  contained  three  data  points  (combinations  of  time  and  dose)  where 
both  injury  level  and  PARP  concentrations  were  measured.  PARP  was  presented  in  cpm/mg  at  2, 
3,  and  5  hours  post-exposure  to  0,  250,  and  500  pM  of  HD.  Level  of  injury  was  presented  as  % 
adherence  at  3,  5,  and  6  hours  following  0,  250,  500,  and  1000  pM  HD  exposures.  In  both  cases, 
we  took  the  ratio  of  PARP  concentration  and  adherence  following  exposure  against  the  control 
values  in  order  to  get  an  estimate  in  the  change  of  PARP  and  adherence  and  then  used  those 
values  as  the  relative  concentration  of  PARP  and  the  relative  injury  level,  respectively.  Plotting 
those  three  points  on  a  graph  and  setting  an  intersection  point  of  zero  (a  value  we  chose  for  the 
baseline  PARP  concentration)  resulted  in  the  curve  shown  in  Figure  6-3. 
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Figure  6-3.  PARP  Concentration  as  a  Function  of  Injury  Level  in  Human  Keratinocytes 
Data  from  (Hinshaw,  Lodhi,  Hurley,  Atkins,  &  Dabrowska,  1999) 


Although  a  quadratic  fit  (using  Excel)  was  the  best  result,  we  chose  a  linear  fit  with  an 
intercept  of  zero  in  order  to  ensure  that:  1)  an  injury  level  of  zero  does  not  result  in  an  undefined 
or  negative  value  2)  PARP  concentrations  start  at  0.0  and  always  rise  as  a  function  of  injury. 
Therefore,  our  equation  for  PARP  concentration  relative  to  injured  cells  is  equal  to: 

PARP  =  5.91  *  NIiX  ( ) 

For  plotting  NAD+  as  a  function  of  PARP,  there  were  six  data  points  that  we  were  able  to 
use,  shown  in  Figure  6-4,  with  the  corresponding  fit.  We  once  again  used  the  relative  PARP 
concentration  values  described  above  and  applied  the  same  calculation  to  the  NAD+ 
concentrations  (presented  in  the  paper  as  pmole/mg)  in  order  to  determine  the  relative  NAD+ 
concentration. 
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Relative  NAD+  Concentrations  as  a  Function  of 
Relative  PARP  Concentrations 
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Relative  PARP  Concentrations 


Figure  6-4.  NAD+  Concentration  as  a  Function  of  PARP  Concentration  in  Human 

Keratinocytes 

We  chose  an  exponential  to  fit  the  NAD/P  ARP  data,  since  the  NAD+  concentration  is 
supposed  to  decrease  as  PARP  increases.  We  needed  to  ensure  that  regardless  of  the  PARP  value 
(which  has  no  bound  on  an  upper  value),  the  NAD+  value  does  not  go  below  zero,  as  that  would 
lead  to  complications  in  the  model.  Fitting  to  an  exponential  function  in  Excel  gives  us  the 
following  equation: 

NAD  =  1.1  *  e~°-51*PARP  (  6-2  ) 

For  plotting  ATP  as  a  function  of  NAD+,  we  only  have  two  data  points  that  we  can  use,  as 
shown  in  Figure  6-5,  with  the  corresponding  fit.  We  used  the  relative  concentration  of  NAD+  as 
described  above  and  applied  the  same  calculation  to  the  ATP  levels  (which  were  presented  in  the 
paper  in  nmole/ 106  cells). 
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Figure  6-5.  ATP  Concentration  as  a  Function  of  NAD+  Concentration  in  Human 

Keratinocytes 


Since  we  only  have  two  useable  data  points,  we  have  no  basis  for  fitting  ATP  as  a  function  of 
NAD+  with  anything  besides  a  straight  line.  Since  NAD+  cannot  go  below  0.0,  and  we  set  the 
intercept  of  the  ATP  function  to  0.0,  the  relative  concentration  of  ATP  will  never  go  below  zero. 
The  fit  gives  us  the  following  equation: 

ATP  =  1.5*  NAD  (6-3) 

Figure  6-6  shows  the  model  estimates  of  PARP,  NAD+,  and  ATP  relative  concentrations  as  a 
function  of  time,  against  the  entire  range  of  possible  injury  levels  (0.0  to  1.0).  The  values  for 
NAD  and  ATP  drop  significantly  as  the  injury  level  increases,  mimicking  the  injury  mechanism 
that  leads  to  increased  necrotic  death,  described  in  the  previous  section. 
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PARP,  NAD+,  and  ATP  Relative  Concentrations 
as  a  Function  of  Injury  Level 
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Injury  Level 

Figure  6-6.  PARP,  NAD+,  and  ATP  as  a  Function  of  Injury 

6.2.2.  Protease  Concentration 

At  the  time  of  writing  this  report,  we  have  not  yet  found  a  study  that  looked  at  both  ATP 
levels  and  protease  activity,  although  we  have  found  several  papers  that  discussed  the  correlation 
between  the  two  (Kehe,  Balszuweit,  Steinritz,  &  Thiermann,  2009)  (Kehe  K.,  2008)  (Pita  & 
Vidal-Asensi,  2010)  (Shakarjian,  et  ah,  2010).  We  do  believe  that  a  more  extensive  literature 
search  can  potentially  lead  to  a  paper  that  explicitly  links  ATP  levels  to  protease  levels,  but  due 
to  time  constraints  on  this  project  we  needed  to  make  use  of  the  literature  we  managed  to  obtain 
and  read  so  far. 

A  recent  paper  (Simbulan-Rosenthal,  et  ah,  2006)  modeled  protease  activity  in  human 
keratinocytes  following  HD  exposure  and  would  have  been  an  excellent  candidate  for  fitting  the 
data  but  the  paper  had  no  doses  and  no  time  points  that  intersected  with  any  of  the  values  used  in 
the  Hinshaw,  et.  al.  study.  A  second  study  (which  we  used  for  the  protease  fit  in  the  HD  model) 
measured  protease  activity  that  was  collected  from  mouse  ear  exposures  (Powers  J.  C.,  1999) 
instead  of  via  in  vitro  human  keratinocyte  exposure.  The  problem  with  using  in  vivo  data  is  that 
we  no  longer  claim  that  the  processes  are  agent-independent,  as  different  agents  will  likely  have 
different  diffusion  rates  in  the  body.  We  instead  need  to  make  the  assumption  that  exposures 
from  all  vesicants  will  result  in  the  same  proteolytic  activity  as  they  would  following  an  HD 
exposure. 
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We  already  analyzed  the  Powers  data  in  the  HD  model  to  determine  the  proteolytic  activity 
increase  over  the  control  as  a  function  of  time.  Table  6-1  shows  the  derived  values  of  proteolytic 
activity  as  a  function  of  time.  This  is  a  fold-increase  over  the  control.  So,  for  example,  elastase 
had  a  0.04  mOD/min  proteolytic  activity  at  hour  three  in  the  absence  of  HD  and  a  0. 102 
mOD/min  proteolytic  activity  after  HD  exposure,  which  results  in  a  2.6-fold  increase  over  the 
control. 


Table  6-1.  Fold  Increase  of  Proteolytic  Activity  Following  HD  Exposure  (derived  from 

Powers  J.  C.,  1999) 


Time 

(hr) 

Elastase 

Tryptase 

Chymase 

Plasmin 

Thrombin 

Cathepsin 

B 

Calpain 

Cathepsin 

H 

3 

2.6 

1.2 

2.9 

4.2 

1.5 

1.3 

1.7 

1.3 

6 

2.1 

1.1 

1.0 

0.2 

1.9 

1.6 

1.7 

1.1 

12 

7.7 

1.3 

0.9 

0.3 

1.3 

1.1 

2.1 

1.2 

24 

19.9 

3.2 

0.5 

4.9 

1.5 

1.0 

2.7 

1.1 

Our  model  does  not  require  us  to  consider  the  activity  of  individual  proteases.  In  the  HD 
model,  we  averaged  these  values  together  to  obtain  one  representative  curve  for  proteolytic 
activity.  Since  we  published  that  paper,  we  have  learned  that  capsase-3  is  considered  to  be  the 
most  important  protease  involved  in  the  mustard  injury  mechanism  (Kehe,  Balszuweit,  Steinritz, 
&  Thiennann,  2009)  (Kehe  K.  ,  2008)  (Simbulan-Rosenthal,  et  ah,  2006).  Described  in  the  table 
above,  caspase-3  activity  was  not  even  measured  in  the  experiment. 

Therefore,  we  had  to  turn  to  the  more  recent  Simbulan-Rosenthal,  et.  al.  (Simbulan- 
Rosenthal,  et  al.,  2006)  study  to  plot  caspase-3  activity  at  16  and  24  hours  as  a  function  of  dose, 
to  first  understand  if  there  is  any  relationship  between  caspase-3  activity  and  dose  and  then  to 
understand  if  there  is  a  relationship  between  caspase-3  activity  and  time.  The  data  from 
Simbulan-Rosenthal  and  the  corresponding  fit  is  shown  in  Figure  6-7  for  16-and  24-hours  post¬ 
exposure.  The  equation  for  caspase-3  activity  as  a  function  of  dose  (D),  as  defined  in  Simbulan- 
Rosenthal  at  16  hours,  is: 


Caspase  —  3  =  0. 024  *  D  —  0.1  (  6~4 ) 

And  the  equation  for  caspase-3  activity,  as  a  function  of  dose  at  24  hours,  is: 

Caspase  —  3  =  0. 03  *  D  +  1  (6-5) 

As  shown  in  the  figure  below,  there  is  strong  linear  relationship  between  caspase-3  activity 
and  dose  at  24-hours  post-exposure  and  a  weaker  relationship  at  16-hours  post-exposure. 
However,  the  best-fit  slopes  are  at  least  of  the  same  order  of  magnitude.  Using  this  equation,  we 
can  now  predict  caspase-3  activity  at  16-  and  24-hours  for  a  250,  500,  and  1000p.M  dose  of  HD, 
the  same  dosage  used  in  the  Hinshaw  study  for  ATP  depletion.  The  result  of  this  calculation  is 
shown  in  Table  6-2. 
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Capsase3  Activity  as  a  Function  of  Dose 


Figure  6-7.  Caspase-3  Activity  as  a  Function  of  HD  Dose  16-  and  24-Hours  Post-Exposure 

At  these  doses,  the  predicted  caspase-3  activity  level  at  16  and  24  hours  are  very  high.  This 
phenomenon  makes  sense,  as  the  injury  mechanism  for  vesicants  states  that,  at  high  doses,  ATP 
depletion  is  more  exaggerated  and  the  level  of  necrosis  increases  significantly.  Also,  the  AMedP- 
8  (C)  injury  profiles  for  HD  show  heavy  blistering  beginning  at  the  18-24  hour  mark.  This  is 
once  again  verified  by  this  model,  since  protease  activity  is  what  leads  to  heavy  blistering.  We 
still  need  to  make  one  more  step;  however,  to  connect  these  values  to  ATP  depletion. 

Table  6-2.  Predicted  Caspase-3  Activity  Levels  for  250,  500,  and  1000  pM  HD  Exposures  at 

16  and  24  Hours  Post-Exposure 


250  pM 

500  pM 

1000  pM 

16  hrs 

5.9 

11.9 

23.9 

24  hrs 

8.5 

16 
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Plotting  the  values  from  the  above  table  and  fitting  the  values  to  a  line  with  an  intercept  of 
1.0  gives  the  plots  shown  in  Figure  6-8. 
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We  can  then  use  the  equations  shown  in  the  figure  above  to  predict  the  caspase-3  activities  at 
5  hours  (the  time  at  which  ATP  levels  were  measured  in  the  Hinshaw  study)  for  250,  500,  and 
1000  pM  exposures.  The  results  of  these  calculations  are  shown  in  Table  6-3,  along  with  the 
relative  ATP  concentrations  derived  from  Hinshaw. 

Table  6-3.  Predicted  Caspase-3  Activity  at  Five  Hours  for  250,  500,  and  1000  pM 

Exposures  of  HD 


250  pM 

500  pM 

1000  pM 

Capsase3  Activity  (5  hours) 

2.55 

4.21 

7.529 

Relative  ATP  Concentration 

0.85 

0.55 

0.5 
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Figure  6-9  shows  the  plot  of  the  table  values  above  and  the  corresponding  fit. 


Figure  6-9.  Caspase-3  Activity  as  a  Function  of  Relative  ATP  Concentration 


Using  the  above  analysis  of  the  Flinshaw  and  Simbulan-Rosenthal  papers,  we  derived  an 
equation  for  protease  that  is  equal  to: 

Protease  =  35.6*  e~3A *ATP  ( ) 

We  realize  that  there  were  several  significant  assumptions  used  to  connect  the  data  from 
these  two  studies.  We  believe  that  a  better  correlation  between  ATP  and  protease  can  be  derived 
if  we  had  a  study  that  explicitly  linked  the  two. 

6.2.3.  Inflammatory  Response 

Once  again,  we  are  lacking  a  comprehensive  study  of  in  vitro  inflammatory  response  as  a 
result  of  killed  cells  and  need  to  derive  the  necessary  data  from  non-ideal  studies.  We  found  one 
study  (Amir,  et  ah,  2000)  that  measured  inflammatory  response  markers  following  FID  exposures 
to  rabbit  eyes.  We  were  able  to  use  the  study  to  make  a  preliminary  correlation  between  killed 
cells  and  inflammatory  response.  The  assumptions  being  made  in  using  this  study  are  that  ocular 
exposures  exhibit  similar  biochemical  processes  to  skin  and  lung  exposures  and  that  rabbit  eyes 
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exhibit  similar  biochemical  processes  to  human  eyes.  For  the  latter  assumption,  the  belief  that 
biochemical  processes  are  system-independent  is  an  underlying  assumption  for  the  entire  model. 
For  the  fonner,  we  do  believe  that  similar  experiments  using  in  vitro  exposures  against  human 
cell  cultures  could  likely  lead  to  different  models.  If  better  data  becomes  available,  we  would 
revisit  the  fit.  However,  our  use  of  the  same  methodology  is  described  below. 

The  Amir  et  al.  study  measures  inflammatory  markers  (protein  levels)  in  the  anterior 
chamber  6-  and  48-hours  post-exposure  of  vapor  HD  against  the  eye.  Additionally,  the  study 
reports  the  size  of  ocular  lesions  (which  can  be  correlated  to  killed  cells)  on  days  1-7  post¬ 
exposure.  The  primary  goal  of  the  study  was  to  measure  the  effects  of  anti-inflammatories 
against  the  ocular  lesions  and  comeal  weight,  but  we  will  only  use  the  data  in  this  section  to 
estimate  the  correlation  between  the  level  of  blistering  (ocular  lesions)  and  the  level  of 
inflammatory  markers  (protein). 

Our  first  step  was  to  fit  the  level  of  corneal  erosions  (reported  in  the  study  as  the  clinical 
score)  as  a  function  of  time,  for  a  dose  of  390|ag/L  over  2  minutes  (the  dose  used  in  the  study). 
The  resulting  plot  and  fit  is  shown  in  Figure  6-10. 
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In  order  to  get  a  good  fit  to  the  data  we  needed  to  use  a  fifth-order  polynomial.  The  equation 
for  ocular  injury  as  a  function  of  time  is  therefore: 

Injury  =  -5E  -  9  *  t5  +  2E  -  6  *  t4  -  0. 0002  *  t3  +  0. 0081  *  t2  +  0. 0087  *  t  +  5E  (  6-7  ) 
-  11 


The  next  step  is  using  Equation  5-7  to  estimate  the  extent  of  corneal  erosions  at  30  minutes,  6 
hours,  and  48  hours;  these  times  are  at  which  the  study  measured  protein  levels  indicative  of 
inflammatory  response.  The  curve  and  resulting  fit  can  be  seen  in  Figure  6-11. 


Figure  6-11.  Inflammatory  Response  as  a  Function  of  Fold  Increase  of  Corneal  Erosions 


The  fit  would  certainly  benefit  from  more  data  points,  but  the  curve  has  the  shape  we  expect 
and  is  a  good  starting  point  for  our  model.  We  see  that  rapid  onset  of  injury  results  in  a 
significant  increase  in  inflammatory  cells.  As  the  injured  cells  decrease,  the  inflammatory 
response  will  decrease  accordingly.  This  would  lead  to  a  rapid  increase  in  inflammatory 
response,  followed  by  a  rapid  drop-off  as  injured  cells  are  vesicated.  Correlating  this  to  the 
AMedP-8  (C)  injury  bands  for  HD,  injury  severity  reaches  a  maximum  in  all  cases  before  48 
hours.  We  infer  that  the  inflammatory  response  should  trend  back  towards  zero  as  there  are  no 
more  injured  and  dead  cells  to  vesicate.  In  order  to  account  for  this  phenomenon,  we  substituted 
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the  exponential  term  with  the  relative  concentration  of  injured  cells.  Our  equation  for 
inflammatory  cells  as  a  function  of  injured  cells  is  therefore: 

Inflam  =.  98  *  exp(Nj2)  (  6-8  ) 

6.2.4.  Calcium  Ion  Concentration 

When  we  developed  the  HD  model  we  stated  that  there  was  no  evidence  that  calcium 
modulation  reduced  the  toxicity  of  HD  injuries.  One  of  the  studies  (Sawer  &  Hamilton,  1999) 
that  supported  this  statement  looked  at  the  effect  that  the  calcium  chelator  BAPTA-AM  and 
calcium  inhibitor  thapsigargin  have  on  cell  viability  following  HD  exposure.  The  study  found 
that  there  were  no  beneficial  therapeutic  effects  from  either  of  these  potential  countermeasures. 
Since  writing  our  final  report  on  the  HD  model,  we  have  come  across  one  paper  that  may  suggest 
a  therapeutic  effect  of  calcium  modulators  (Ray,  Benton,  Anderson,  Byers,  &  Petrali,  2000), 
although  not  in  a  way  that  is  consistent  with  our  model.  The  findings  of  Ray  et  al.  suggest  that 
calcium  modulators  may  actually  work  by  slowing  cellular  proliferation  and  metabolic  rates. 
Thus,  the  rate  of  apoptotic  cell  death  slows,  but  not  necrotic  cell  death,  as  we  have  stated  in  our 
model. 

There  are  the  remaining  questions  of  whether  or  not  the  dosing  of  calcium  modulators 
required  for  a  therapeutic  effect  is  toxic  in  an  organism,  or  whether  or  not  slowing  the  rate  of 
apoptotic  cellular  death  would  actually  stop  the  formation  of  blisters  that  follow  large  exposures 
of  vesicants.  The  author  references  having  begun  experimentation  in  animals  using  BAPTA-AM, 
but  we  have  not  found  a  follow-up  report  with  the  results. 

Pending  future  data,  we  still  consider  calcium  ions  and  calcium  modulation  to  not  be  integral 
components  of  vesicant  injury.  We  have  defined  a  simple  equation  for  calcium  ions  that  is 
independent  of  injury  and  dose  and  will  have  no  effect  on  the  model  output: 

Ca  =  1  ( 6- 9) 

In  our  model,  calcium  ions  will  be  held  at  a  constant  concentration.  If  further  research  begins 
to  show  a  better  understanding  of  calcium  ions’  role  in  vesication  or  the  therapeutic  effects  of 
calcium  ions,  we  will  revisit  this  assumption. 

6.2.5.  Glutathione  (GSH)  Suppression 

We  found  two  excellent  papers  that  studied  the  effects  several  vesicants  have  on  glutathione 
levels  in  various  compartments  of  the  body. 

The  first  study  (Sharma,  Vijayaraghavan,  &  Agrawal,  2010)  looked  at  the  effects  HD,  HN-1, 
HN-2,  and  HN-3  have  on  body  weight,  spleen  weight,  white  blood  cells,  red  blood  cells,  liver 
enzymes,  GSH,  oxidized  glutathione  (GSSG),  liquid  peroxidation  (MDA),  and  DNA  damage. 
The  effects  are  followed  in  percutaneous  exposure  to  mice  ears  with  and  without  treatment.  This 
is  a  model  study  in  the  type  of  data  that  can  best  be  used  for  the  GVM,  because  it  considers 
different  vesicants,  different  compartments,  and  different  countermeasures.  Additionally,  by 
focusing  on  systemic  effects,  it  allows  us  to  model  GSH  suppression  via  systemic  dose.  The  first 
point  to  make  after  reviewing  this  paper  is  that  HN-1,  HN-2,  HN-3,  and  HD  all  result  in 
approximately  a  60%  decrease  in  GSH  levels  in  the  liver.  This  fact  further  strengthens  two  of  our 
assumptions:  the  underlying  biological  process  disruptions  are  largely  agent-dependent  and  our 


37 


assumption,  the  three  nitrogen  mustard  analogues  are  as  toxic  as  sulfur  mustard.  In  this  study, 
GSH  suppression  was  only  measured  in  the  liver,  but  a  second  study  (Lomash,  Deb,  Rai,  Jadhav, 
R.Vijayaraghavan,  &  S.C.Pant,  2011)  gives  some  insight  into  why. 

The  Lomash  et  al.  study  examined  the  effect  of  HD  on  body  weight,  liver,  kidney,  and  spleen 
weight,  white  blood  cells,  red  blood  cells,  hemoglobin,  GSH  levels  in  the  blood  and  liver,  GSSG 
and  MDA  levels  in  the  liver,  and  even  the  occurrence  and  severity  of  injury  in  the  liver,  kidneys, 
spleen,  and  skin.  The  purpose  of  this  extensive  study  was  to  design  an  in  vivo  model  that 
produced  skin  lesions  without  producing  systemic  toxicity.  As  a  result,  the  study  gave  significant 
insight  into  the  role  of  GSH  suppression  in  vesicant  injury.  This  paper  benefits  our  model  by 
providing  evidence  that  GSH  plays  a  significant  role  in  systemic  toxicity.  Small  exposures  of  HD 
(5mg/kg)  result  in  little  to  no  systemic  toxicity,  but  still  produce  significant  skin  lesions,  whereas 
larger  exposures  of  HD  (lOmg/kg)  result  in  systemic  toxicity  in  almost  all  animals  exposed  (and 
also  produces  significant  skin  lesions).  The  second  insight  gleaned  from  this  study  is  that  GSH 
levels  in  the  blood  are  not  significantly  changed  when  challenged  with  a  5  mg/kg  or  lOmg/kg 
percutaneous  dose.  However,  GSH  levels  in  the  liver  are  significantly  changed  when  challenged 
with  a  10  mg/kg  dose,  but  not  with  a  5  mg/kg  dose.  This  leads  us  to  believe  that  it  is  GSH- 
suppression  in  the  target  organs  that  results  in  systemic  toxicity  and  not  in  the  blood.  This  study 
would  have  helped  further  if  the  researchers  measured  GSH  suppression  in  the  bone  marrow  and 
spleen,  as  we  would  potentially  be  able  to  create  a  toxicodynamic/kinetic  model  of  vesicant 
distribution  and  effect.  For  now  we  have  to  group  all  of  the  important  organs  into  one 
compartment  and  apply  the  systemic  dose  to  that  compartment.  The  only  data  that  we  can 
correlate  dose-to-effect  is  in  the  liver.  Since  this  is  a  new  research  area  (the  two  papers  were 
published  in  2010  and  2011,  respectively),  we  are  hopeful  that  continued  research  allows  us  to 
better  refine  this  model. 

Another  benefit  of  the  study,  is  that  we  can  extrapolate  the  findings  to  human  exposures 
using  the  toxicokinetic  model  we  are  implementing  in  the  GVM. 

First  we  need  to  calculate  the  systemic  dose  that  results  from  a  5  mg/kg  and  a  lOmg/kg 
percutaneous  exposure  in  a  70-kg  man  (350  and  700mg,  respectively).  Using  that  as  input  to  the 
toxicokinetic  algorithm,  we  get  the  two  estimates  of  systemic  dose  as  a  function  time,  shown  in 
Figure  6-12. 


38 


Systemic  Dose  Over  Time  for  Two  Percutaneous 


Figure  6-12.  Systemic  Dose  Over  Time  for  a  350mg  and  700mg  Percutaneous  Exposure 


The  Sharma  et  al.  study  measured  GSH  levels  in  the  liver  on  days  4  and  8.  As  shown  in  the 
figure  above,  most  of  the  systemic  dose  enters  the  circulation  system  before  day  2.  The  Lomash 
et  al.  study  measures  GSH  levels  in  the  liver  and  kidneys  on  days  3  and  7.  Table  6-4  tabulates  the 
data  from  the  two  studies  (Sharma  et  al.)  using  a  percutaneous  dose  of  8. 1  mg/kg. 

Table  6-4.  GSH  Supression  in  the  Liver  for  Three  Different  Doses 


Percutaneous  dose 

5  mg/kg 

8.1  mg/kg 

10  mg/kg 

3  days 

1 

- 

0.8 

4  days 

- 

0.4 

- 

7  days 

1 

- 

0.8 

8  days 

- 

0.4 

- 

One  issue  is  immediately  apparent:  Sharma  et  al.  study  measured  a  greater  level  of  GSH 
suppression  in  the  liver  for  an  8.1  mg/kg  dose  than  the  Lomash  et  al.  study  measured,  using  a  10 
mg/kg  dose.  This  is  likely  due  to  one  big  difference  between  the  two  studies:  For  the  Shanna  et 
al.  study  the  mice  were  exposed  on  their  backs;  however,  for  the  Lomash  et  al.  study  the  mice 
were  exposed  on  their  ears.  This  tends  to  be  an  issue  when  working  with  percutaneous  hazards; 
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effects  tend  to  be  exposure-site  dependent.  In  this  case  we  can  postulate  several  reasons  for  the 
drastic  difference  between  the  two  exposure  sites: 

1 .  Differences  in  thickness  between  back  tissue  and  ear  tissue 

2.  Differences  in  fat  content  between  the  back  tissue  and  the  ear 

3.  Differences  in  the  paths  required  to  get  to  the  target  tissue,  especially  if  there  are  multiple 
sinks  or  metabolites  en  route 

Regardless  of  the  reason,  the  issue  highlights  one  of  the  problems  with  modeling 
percutaneous  exposures.  For  this  implementation,  we  assume  that  a  human  exposure  would 
resemble  the  ear  exposure  (Lomash  et  ah),  considering  we  are  modeling  to  a  70-kg  warfighter 
who  would  most  likely  not  have  an  exposed  back.  From  a  mathematical  standpoint,  the  Lomash 
study  also  gives  us  two  data  points  to  include  in  the  model.  The  plot  and  corresponding  fit  are 
shown  in  Figure  6-13.  The  equation  for  GSH  suppression  as  a  function  of  systemic  dose  is: 

GSH  =  1.25  *  e-°OO7*0s  (6-10) 

We  chose  to  fit  GSH  using  an  exponential  curve,  so  that  the  LD50  for  percutaneous 
exposures  results  in  a  positive  value  for  GSH  suppression.  For  example,  using  our  model,  a  1400 
mg  exposure  of  HD  (the  LD50  for  percutaneous  exposures)  produces  a  126mg  systemic  dose  and 
a  GSH  level  of  0.5 1 .  We  will  discuss  how  we  tie  GSH  suppression  to  lethality  in  the  coming 
sections. 
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6.3  The  Effect  of  Biological  Processes 

Now  that  we  have  defined  the  models  that  act  as  the  foundation  of  the  GVM,  we  can  overlay 
the  effects  that  perturbing  these  processes  have  on  cells.  There  are  four  primary  categories  for 
cells  in  our  model:  normal  healthy  cells,  injured  cells,  necrotic  cells,  and  apoptotic  cells.  We  will 
first  define  the  equations  and  then  complete  the  fits. 

6.3.1.  Proportion  of  Healthy  Cells  (Nn) 

The  proportion  of  healthy  cells  in  the  basal  epidermis  and  dermis  is  given  by  the  variable  Nn. 
It  describes  the  proportion  of  cells  at  an  injury  site  that  have  undamaged  DNA  with  respect  to 
vesicant  injury. 

At  equilibrium  (pre-exposure),  we  assume  the  value  of  Nn  is  equal  to  1.0,  meaning  all  of  the 
cells  at  the  injury  site  have  not  been  damaged  by  vesicant  exposure.  Nn  is  modified  by: 

1 .  The  rate  of  mitotic  cell  division,  which  we  assume  to  be  proportional  to  the  rate  of 

programmed  cell  death.  Therefore,  pre-exposure,  the  rate  of  mitotic  division  results  in  no  net 
gain  of  healthy  cells.  Post-exposure,  the  rate  of  mitotic  cell  division  will  result  in  an  increase 
in  the  value  of  Nn,  until  equilibrium  is  once  again  reached. 
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2.  The  rate  of  DNA  alkalization  from  vesicant  exposure.  A  cell  with  damaged  DNA  is  moved 
into  the  cohort  Ni,  which  represents  the  proportion  of  cells  with  injured  DNA. 

3.  The  rate  of  DNA  repair.  A  cell  with  damaged  DNA  can  once  again  move  back  to  the  Nn 
cohort,  provided  there  are  significant  levels  of  PARP  and  NAD+  to  repair  the  DNA. 
Although  research  has  shown  that  PARP  and  NAD+  can  incorrectly  repair  DNA,  we  are  not 
modeling  that  mechanism  or  its  effect  in  this  implementation  of  the  model. 

We  can  now  draw  a  system  diagram  for  percutaneous  injuries  that  accounts  for  liquid  and 
vapor  vesicant,  as  shown  in  Figure  6-14. 


Figure  6-14.  Proportion  of  Healthy  Cells  in  the  Skin  (Nn) 


The  equation  that  describes  the  above  diagram  is: 

d  (6-H) 

—  —  ksp;it  *  (1  Nn)  *  Nn  kbind  *  D  *  Nn  +  krepair  *  PARP  *  NAD  *  Nj 

Where: 

1 .  D  is  the  dose  in  the  skin,  given  in  milligrams  and  described  in  Section  5, 

2.  kgpiit  is  the  rate  of  mitotic  division  given  in  min'1, 

3.  kbmd  is  the  rate  vesicant  binds  with  the  DNA  in  healthy,  normal  cells  in  min'1, 

4.  krepair  is  the  rate  of  DNA  repair  given  in  min'1,  and 

5.  NAD  and  PARP  are  the  relative  concentrations  of  NAD+  and  PARP  in  the  skin,  respectively, 
and  are  defined  by  the  equations  in  Section  5.2.1. 

6.3.2.  Proportion  of  Injured  Cells  (Ni) 

The  proportion  of  injured  cells  is  given  by  the  variable  Ni.  It  describes  the  proportion  of  cells 
at  an  injury  site  that  have  fragmented  DNA  with  respect  to  vesicant  injury. 
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At  equilibrium  (pre-exposure),  we  assume  the  value  of  Ni  is  equal  to  0.0,  meaning  none  of 
the  cells  at  the  injury  site  have  been  damaged  by  vesicant  exposure.  Ni  is  modified  by: 

1 .  The  rate  of  mitotic  cell  division,  which  we  assume  to  be  proportional  to  the  rate  of 
programmed  cell  death.  For  injured  cells,  mitotic  division  causes  the  cell  to  die.  For  the 
percutaneous  model,  this  specifically  places  the  cell  in  the  apoptotic  compartment. 

2.  The  rate  of  DNA  alkalization  from  vesicant  exposure.  A  cell  with  damaged  DNA  is  moved 
into  the  injured  cohort. 

3.  The  rate  of  DNA  repair.  An  injured  cell  can  once  again  move  back  to  the  Nn  cohort, 
provided  there  are  significant  levels  of  P ARP  and  NAD+  to  repair  the  DNA. 

4.  The  rate  of  necrotic  cell  death.  For  the  inhalation  and  ocular  model,  the  principle  behind  this 
rate  in  the  GVM  is  the  same  as  the  HD  model;  an  injured  cell  dies  at  a  constant  rate,  provided 
it  wasn’t  repaired  by  PARP  and  NAD+.  For  the  percutaneous  model,  this  rate  of  necrotic 
cellular  death  is  tied  to  the  depletion  of  ATP,  which  increases  the  concentration  of  protease, 
as  described  in  section  5.2 

We  can  now  draw  a  system  diagram  for  percutaneous  injuries  that  accounts  for  liquid  and 
vapor  vesicant,  as  shown  in  Figure  6-15. 
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Figure  6-15.  Proportion  of  Injured  Cells  in  the  Skin  (Nfi 


The  equation  that  describes  the  above  diagram  is: 

d  (6-12) 

—  NI  —  kbind*D  *  Nn  kspUt  *  (1  Nj)  *  Nj  krepair  *  PARP  *  NAD  *  NI  kkill 

*  Prot  *  Ca*  Inflam  *  N} 


Where: 

1 .  D  is  the  percutaneous  dose,  given  in  milligrams, 

2.  kSpiit  is  the  rate  of  mitotic  division  given  in  min’1, 
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3.  kbmd  is  the  rate  DNA  vesicant  with  the  DNA  in  healthy,  normal  cells  in  min'1, 

4.  krepair  is  the  rate  of  DNA  repair  given  in  min'1, 

5.  kkin  is  the  rate  of  DNA  death  given  in  min'1, 

6.  NAD,  PARP,  Prot,  Ca,  and  Inflam  are  the  relative  concentrations  of  NAD+,  PARP,  protease, 
calcium  ions,  and  inflammatory  cells,  respectively,  as  described  in  Section  5.2 

6.3.3.  Proportion  of  Apoptotic  Cells  (Nk,a) 

The  proportion  of  cells  undergoing  apoptotic  cell  death  is  given  by  the  variable  Nk,a-  It 
describes  the  proportion  of  cells  at  an  injury  site  that  die  as  a  result  of  damaged  DNA  attempting 
mitotic  cell  division.  We  only  consider  apoptotic  cell  death  for  percutaneous  injuries  in  this 
model.  The  only  data  we’ve  come  across  that  could  support  an  apoptotic/necrotic  model  and  give 
weight  to  the  benefit  of  differentiating  has  been  done  in  skin  cultures  and  via  percutaneous 
exposures  on  animals.  It  is  possible  we  may  have  enough  data  to  justify  an  apoptotic/necrotic 
model  for  the  eyes  and  lung,  at  some  point,  but  have  yet  to  find  it  at  the  time  of  writing  this 
report. 

At  equilibrium  (pre-exposure),  we  assume  the  value  of  Nk,a  is  equal  to  0.0,  meaning  none  of 
the  cells  at  the  injury  site  have  been  killed  via  apoptotic  cellular  death,  resulting  from  vesicant 
exposure.  This  is  an  important  distinction,  as  the  primary  mode  of  cellular  death  is  via  the 
apoptotic  pathway.  The  only  factor  that  contributes  to  apoptotic  cell  death  in  our  model  is: 

1 .  The  rate  of  mitotic  cell  division,  which  causes  injured  cells  to  die  in  our  model.  For  the 
percutaneous  model,  this  specifically  places  the  cell  in  the  apoptotic  compartment. 

We  can  draw  a  system  diagram  for  apoptotic  cell  death,  as  shown  in  Figure  6-16. 
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Figure  6-16.  Proportion  of  Apoptotic  Cells  in  the  Skin  (Nk,a) 


The  equation  that  describes  the  above  diagram  is: 

d 

—  kspnt  *  (1  —  Nj)  *  Nj 

Where: 

1 .  kspiit  is  the  rate  of  mitotic  division  given  in  min'1,  and 

2.  Ni  is  the  proportion  of  injured  skin  cells 


(6-13) 
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6.3.4.  Proportion  of  Necrotic  Cells  (Nk,n) 

The  proportion  of  cells  undergoing  necrotic  cell  death  is  given  by  the  variable  Nk,n-  It 
describes  the  proportion  of  cells  at  an  injury  site  that  die  as  a  result  of  damaged  DNA  being 
repaired  by  PARP.  Thus,  NAD+  is  depleted,  which,  in  turn,  depletes  ATP  stores. 

At  equilibrium  (pre-exposure),  we  assume  the  value  of  Nk,n  is  equal  to  0.0,  meaning  none  of 
the  cells  at  the  injury  site  have  been  killed  via  necrotic  cellular  death.  The  factors  that  contribute 
to  necrotic  cell  death  in  our  model  are: 

1.  The  rate  of  mitotic  cell  division  in  injured  ocular  and  lung  cells.  Attempted  mitotic  division 
in  injured  skin  cells  results  in  apoptotic  cell  death. 

2.  The  rate  of  cell  death  due  to  ATP  depletion. 

We  can  now  draw  a  system  diagram  for  necrotic  cell  death  in  the  skin,  as  shown  in  Figure 
6-17. 


Figure  6-17.  Proportion  of  Necrotic  Cells  in  the  Skin  (Nn,k) 


The  equation  that  describes  the  above  diagram  is: 

d  (6-14) 

-j-^  Nr.n  =  kkiii  *  Pfot  *  Inf  lam  *  Ca*  N{ 

Where,  k^iii  is  the  rate  of  necrotic  cell  death  given  in  min'1 . 

6.4  Fitting  Cell  Concentration  Equations 

We  are  getting  to  the  point  in  the  model  where  we  will  need  to  start  differentiating  between 
the  different  agents,  primarily  for  the  value  kbind,  which  defines  the  rate  at  which  vesicant  binds 
with  and  fragments  DNA.  This  is  an  agent-dependent  variable.  However,  we  are  still  making  the 
assumption  that  the  other  three  rates,  kspiit,  krepair,  and  kkin,  are  still  agent-independent  for  reasons 
specified  below: 

•  kspiit  is  the  rate  of  mitotic  cell  division,  a  biological  process  that  changes  with  the  level  of 
injury,  regardless  of  how  the  injury  was  caused. 
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•  krepair  is  the  rate  at  which  PARP  repairs  fractured  DNA  and  is  independent  of  how  the  DNA 
became  injured. 

•  kyii  is  the  rate  at  which  protease  breaks  down  injured  cells  in  the  absence  of  ATP  and  is 
independent  of  how  the  ATP  was  depleted  or  the  cell  was  injured. 

Due  to  this  assumption,  we  fit  all  four  parameters  using  HD  as  a  baseline  and  then  made  the 
necessary  adjustments  to  kbind  for  the  other  agents. 

6.4.1.  The  Rate  of  Mitotic  Division  (kspiit) 

We  discussed  the  rate  of  mitotic  division  in  the  HD  paper.  For  that  model,  we  let  the  rate  of 
cellular  division  ultimately  be  a  fit  parameter.  However,  for  the  GVM,  we  set  the  rate  of  mitotic 
division  at  0.004hf 1  in  order  to  match  the  value  in  the  literature  (Murray,  Mathematical  Biology, 
Volume  2,  1993).  We  assume  that  this  rate  is  the  same  for  all  three  exposure  sites  of  interest.  We 
are  assuming  that  this  rate  is  constant  for  all  exposure  sites. 

6.4.2.  The  Repair  Rate,  Kill  Rate,  and  HD  Binding  Rate  for  Percutaneous  Exposures 

After  setting  the  rate  of  mitotic  division  to  a  rate  found  in  the  literature  we  are  left  with  three 
equations  and  three  unknowns,  an  easy  problem  to  solve  using  an  equation  solver  like  Berkeley 
Madonna.  As  mentioned,  we  first  calibrated  the  model  to  HD,  as  this  would  give  us  values  for 
the  repair  rate  and  kill  rate  in  each  of  the  four  systems.  We  can  then  calibrate  the  binding  rate  to 
each  of  the  remaining  agents.  For  all  fits  we  are  trying  to  capture  the  rate  of  injury  and  repair 
documented  in  AMedP-8  (C).  For  the  HD  model  we  attempted  to  match  AMedP-8  (C)  severity 
scales  to  the  proportion  of  injured  and  killed  cells,  but  for  the  GVM  we  assume  that  is  an 
unnecessary  step.  The  AMedP-8  (C)  severity  scales  have  no  real  correlation  to  magnitude  of 
injury  (i.e.,  a  Severity  Scale  of  2  is  not  necessarily  “twice  as  a  severe”  as  a  Severity  Scale  of  1). 
Mapping  a  quantitative  value  to  a  qualitative  description  of  illness  is  not  a  valid  mechanism. 

For  the  GVM  we  define  an  Injury  function,  which  is  the  complement  of  the  normal,  healthy 
cell  concentration: 


Injury  =  1  —  NNx  (  6-15 ) 

In  doing  so,  we  are  stating  that  the  level  of  injury  is  directly  related  to  the  concentration  of 
nonnal,  healthy  cells.  If  an  exposure  site  has  a  1 .0  concentration  of  nonnal,  healthy  cells  the 
model  predicts  that  there  is  no  injury.  If  an  exposure  site  had  a  0.5  concentration  of  nonnal, 
healthy  cells  the  model  predicts  that  only  half  of  the  exposure  site  is  operational.  We  fit  the 
Injury  equation,  so  that  it  has  a  similar  peak  and  duration  as  AMedP-8  (C)  severity  scales. 

6.4.3.  Skin  Injury 

In  order  to  give  an  example  of  the  fitting  process,  we  focus  on  percutaneous  exposures.  The 
top  dose  band  for  percutaneous  exposures  in  AMedP-8  (C)  is  for  exposures  greater  than  180  mg- 
min/m3-  As  shown  in  Figure  6-18  (green  line),  an  individual  reaches  Severity  Level  1  at  2  hours, 
Severity  Level  2  at  5  hours,  and  Severity  Level  3  at  24  hours.  Severity  Level  3  is  the  injury  peak. 
On  the  recovery  side,  an  individual  goes  back  down  to  a  Severity  Level  1  at  504  hours  and  is 
fully  recovered  at  588  hours.  Constraining  the  equation  for  Injury  to  meet  these  criteria,  we  got 
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the  parameters  for  percutaneous  exposure  to  HD  that  are  shown  in  Figure  6-18.  These  values 
were  derived  using  the  curve  fitting  function  in  Berkeley  Madonna. 


AMedP-8  (C)  HD  S/S  Severity  -  Percutaneous  Exposure 


Figure  6-18.  AMedP-8  (C)  HD  S/S/  Severity  for  Percutaneous  Exposures 


Table  6-5.  HD  Percutaneous  Values  for  kbmd,  krepair,  and  kkni 


kbind  (hr  ) 

^repair  (hr  ) 

kkiu  (hr'1) 

0.04 

0.003 

0.003 
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Using  the  values  in  the  table  above,  we  get  the  results  shown  in  Figure  6-19. 


Figure  6-19.  Percutaneous  Model  Output  for  HD  Compared  to  AMedP-8  (C) 


Note  how  the  model  peaks  at  the  exact  time  as  the  injury  profile.  Recovery  is  also  consistent 
with  the  AMedP-8  (C)  output.  For  the  severe  percutaneous  dosage  we  get  a  maximum  value  of 
0.7  for  the  Injury  function,  meaning  70%  of  the  cells  in  the  dermis  and  basal  epidermis  are 
injured. 

Figure  6-20  shows  the  concentrations  of  ATP,  PARP,  NAD+,  and  protease,  following  a  180 
mg/m  exposure.  Note  the  rapid  increase  in  PARP  that  leads  to  a  rapid  decrease  in  NAD+  and 
ATP,  which  causes  a  rapid  increase  in  protease. 
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Underlying  Biological  Processes  for  a 


0.1  1  10  100  1000 

Time  (Hours) 

Figure  6-20.  ATP,  PARP,  NAD+,  and  Protease  Concentrations  for  a  180mg/m2 

Percutaneous  Exposure  to  HD. 

2 

Figure  6-21  shows  the  percentage  of  necrotic  and  apoptotic  cells,  following  a  180mg/m 
exposure  of  HD.  Note  that  at  these  high  exposures,  the  majority  of  injured  cells  become  necrotic, 
indicating  a  severe  injury.  The  model  equations  were  designed  in  such  a  way  that  would  make 
apoptotic  cell  death  the  primary  mechanism  for  low  exposures  and  necrotic  cell  death  the 
primary  mechanism  for  high  exposures.  Figure  6-22  shows  this  phenomenon.  A  1  Omg/rrT 
exposure  leads  to  a  higher  percentage  of  apoptotic  cell  death.  As  the  exposure  increases,  the 
apoptotic  pathway  begins  to  apex,  while  the  necrotic  pathway  becomes  the  dominant  pathway. 
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Figure  6-21.  Relative  Concentration  of  Necrotic  vs.  Apoptotic  Cell  Death  for  a  180mg/m2 

Percutaneous  Exposure  to  HD. 
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Figure  6-22.  Relative  Concentration  of  Necrotic  vs.  Apoptotic  Cell  Death  for  Four  Different 

Percutaneous  Exposures  of  HD. 
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6.4.4.  HN-1,  HN-2,  HN-3,  and  Lewisite  Binding  Rate 

The  difficulty  in  modifying  the  binding  rates  for  these  four  agents  is  the  fact  that  FM  3-1 1.9 
reports  that  they  have  the  same  toxicities  as  HD.  This  leaves  us  with  two  choices: 

1 .  Assume  that  all  five  agents  have  the  same  exact  binding  rates. 

2.  Use  the  published  EPA  AEGL  values  (which  are  different  for  some  of  the  agents)  to  scale 
these  binding  rates. 

We  decided  to  go  with  the  second  option,  since  the  AEGL  values  represent  a  more  recent 
review  of  agent  toxicity.  The  AEGL  values  for  the  five  vesicants  are  shown  in  Figure  6-23. 


AEGL-3  Values  for  Five  Different  Vesicants 


Figure  6-23.  AEGL-3  Concentrations  for  HD,  HN-1,  HN-2,  HN-3,  and  Lewisite 


For  the  purpose  of  this  model,  we  are  most  concerned  with  the  10-minute  AEGL-3 
concentrations.  Lewisite  and  HD  have  equivalent  10-minute  AEGL-3  values,  adding  support  to  a 
decision  to  assume  that  lewisite  and  HD  have  equivalent  binding  rates.  HD  and  lewisite  do 
deviate  from  each  other  with  longer  exposures,  but  we  do  not  yet  have  a  mechanism  to  account 
for  this  deviation.  It  is  possible  that  the  difference  could  come  from  the  variances  in  the 
toxicokinetics  of  the  two  agents;  perhaps,  lewisite  metabolizes  more  rapidly  than  HD.  Due  to 
time  and  data  constraints,  we  were  not  able  to  account  for  this  difference  in  this  implementation 
of  the  model. 
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•3 

The  AEGL  documentation  does  state  that  being  exposed  to  a  2.2  mg/m  concentration  of 
nitrogen  mustard  for  ten  minutes  will  have  an  equivalent  effect  as  being  exposed  to  a  3.9  mg/m3 
concentration  of  HD  for  ten  minutes,  making  nitrogen  mustard  more  toxic  than  HD.  We  do  not 
believe  this  difference  can  be  ignored.  In  order  to  account  for  this  difference  using  the  binding 
rate,  we  ran  the  inhalation  model  using  a  3.9  mg/m  exposure  for  ten  minutes.  This  resulted  in  a 
peak  injury  level  of  approximately  0.3 1 .  We  wanted  to  derive  a  binding  rate  for  nitrogen  mustard 
that  resulted  in  the  same  injury  level  as  a  2.2  mg/m  exposure.  This  was  achieved  with  a  binding 
rate  of  0.6  hr"1.  Since  this  binding  rate  represents  an  increase  of  1.5  for  inhalation  exposures,  we 
made  the  assumption  that  the  binding  rate  for  the  other  three  injury  sites  are  also  1.5  times 
greater.  All  of  the  rate  parameter  values  are  summarized  in  Table  6-6. 

6.5  Summary  of  Rate  Parameters 


Table  6-6.  Summary  of  all  Rate  Parameter  Values 


Ocular 

Inhalation 

Percutaneous 

ksplit 

0.004 

0.004 

0.004 

kbind?  HD 

0.4 

0.35 

0.04 

kbind,  HN-1,  -2,  and  -3 

0.6 

0.6 

0.06 

kbind,  lewisite 

0.4 

0.4 

0.04 

Iz 

^repair 

10 

10 

0.003 

k-kill 

5.8 

2 

0.003 

6.6  Calculating  Casualties 

Since  the  vesicant  model  was  calibrated  using  AMedP-8  (C)  severity  profdes,  we  need  to 
make  use  of  the  dose  bands  presented  in  the  reference  document  when  calculating  casualties.  We 
correlate  the  dose  bands,  defined  in  AMedP-8  (C)  to  the  proportion  of  killed  cells  (necrotic  cells, 
particularly,  for  percutaneous  exposures),  to  determine  the  severity  level  for  a  given  exposure. 
For  example,  Table  6-7  shows  the  range  of  AMedP-8  (C)  dose  bands  and  the  corresponding 
proportion  of  necrotic  cells  that  are  calculated  in  the  model.  To  calculate  an  individual’s  severity 
profile  from  AMedP-8  (C),  you  would  calculate  the  proportion  of  necrotic  cells  for  a  given 
dosage  and  see  where  the  values  lies  in  relation  to  Table  6-7. 
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Table  6-7.  AMedP-8  (C)  Percutaneous  Dose  Ranges  Compared  to  Proportion  of  Necrotic 
Cell  Ranges  Following  Percutaneous  Exposure 


AMedP-8  (C)  Dose  Range 

Proportion  of  Necrotic  Cells  Range 

12  -  <  125  mg-min/m3 

0.01  -<0.37 

125  -  <  180  mg-min/m3 

0.37 -<0.5 

>180  mg-min/m3 

>0.5 

3 

So  if  an  individual  is  exposed  to  100  mg-min/m  of  HD,  he  would  have  a  proportion  of 
necrotic  cells  equal  to  0.39.  This  value  falls  in  between  0.37  and  0.5  and  therefore,  the 
individual’s  sign/symptom  severity  profile  follows  the  125  -  <180  mg-min/m  profile  from 
AMedP-8  (C).  As  we’ll  show  in  the  next  section,  certain  countermeasures  will  reduce  the 
proportion  of  necrotic  cells  and  could  potentially  move  an  individual  into  a  lower  AMedP-8  (C) 
dose  band. 

Using  the  same  method  detailed  above,  we  calculated  the  proportion  of  killed  cells  dose 
range  for  ocular  (Table  6-8)  and  inhalation  (Table  6-9)  exposures. 

Table  6-8.  AMedP-8  (C)  Ocular  Dose  Ranges  Compared  to  Proportion  of  Killed  Cell 

Ranges  Following  Ocular  Exposure 


AMedP-8  (C)  Dose  Range 

Proportion  of  Killed  Cells  Range 

4  -  <  26  mg-min/m3 

0.05  -<0.21 

26  -  <  50  mg-min/m3 

0.21  -<0.38 

50  -  <  70  mg-min/m3 

0.38  -<.51 

70  -  <  100  mg-min/m3 

0.51  -<0.73 

>100  mg-min/m3 

>0.73 

Table  6-9.  AMedP-8  (C)  Inhalation  Dose  Ranges  Compared  to  Proportion  of  Killed  Cell 

Ranges  Following  Inhalation  Exposure 


AMedP-8  (C)  Dose  Range 

Proportion  of  Necrotic  Cells  Range 

50  -  <  70  mg-min/m3 

0.16  -<0.22 

70  -  <  100  mg-min/m3 

0.22  -<0.3 

100  -  <  150  mg-min/m3 

0.3  -<0.45 

150  -  <  250  mg-min/m3 

0.45  -<0.8 

250  -  <  1200  mg-min/m3 

0.8  -<  1.0 

>  1200  mg-min/m3 

>  1.0 
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6.7  Calculating  Fatalities 

In  AMedP-8  (C),  a  fatality  is  someone  who  spends  15  minutes  at  Severity  Level  4.  This 
would  make  everyone  who  receives  an  inhalation  dose  greater  than  250  mg-min/m  ,  a  fatality. 
Using  the  criteria  for  casualty  calculation,  we  can  adopt  the  same  methodology  in  the  GVM  for 
inhalation  exposures.  If  an  individual  spends  greater  than  10  minutes  with  a  proportion  of 
necrotic  cells  greater  than  0.8,  the  individual  becomes  a  fatality. 

For  systemic  effects,  however,  we  must  make  use  of  the  GSH  suppression  level,  since  an 
LD50  for  percutaneous  exposures  (1400  mg/man)  will  lead  to  complete  necrosis  of  the  injury  site 
as  would  significantly  lower  doses.  Therefore,  necrosis  at  the  injury  site  is  not  a  good  metric  for 
systemic  lethality  in  the  GVM.  AMedP-8  (C)  uses  >  1400  mg  per  man  as  the  fatal  percutaneous 
band.  This  value  is  equivalent  to  the  LD50  in  FM  3-1 1.9.  Since  we  don’t  have  any  other  systemic 
dosage  bands  in  AMedP-8  (C),  we  decided  to  make  use  of  the  toxicity  data  for  percutaneous 
exposures  for  calculating  the  incidence  of  lethality  via  systemic  injuries. 

FM  3-1 1.9  specifies  an  LD5o  of  1400  mg  per  man  and  a  probit  slope  of  7  for  HD,  FIN-1,  HN- 
2,  HN-3,  and  lewisite.  Using  these  values  to  calculate  20  doses  between  the  LD01  (651  mg/man) 
and  the  LD99  (3010  mg/man)  and  using  those  doses  as  inputs  into  the  GVM  results  in  a  range  of 
GSH  suppression  levels  from  0.17  to  0.81.  Figure  6-24  shows  the  plot  of  incidence  of  fatality  as 
a  function  of  GSH  suppression,  along  with  the  corresponding  fit  to  a  log-normal  distribution 
(dashed  line). 


Incidence  of  Fatality  as  a  Function  of  GSH 
Concentration 


0  0.2  0.4  0.6  0.8  1  1.2 

GSH  Concentration 


Figure  6-24.  Incidence  of  Fatality  as  a  Function  of  GSH  Concentration  Following  Systemic 

Exposures. 


54 


The  equation  for  the  log-normal  distribution  of  Probability  of  Fatality  for  a  given  GSH 
suppression  level  is: 

1  —  ln(l  —  GSH)  +  p  (6-16) 

Probability  of  Fatality  =  -erfc( - — - ) 

2  a  *  V2 


Where: 

1 .  erfc  is  the  complementary  error  function, 

2.  GSH  is  the  concentration  of  GSH, 

3.  p  is  the  mean,  equal  to  0.74,  and 

4.  a  is  the  standard  deviation,  equal  to  0.29. 

In  order  to  calculate  whether  or  not  someone  becomes  a  fatality,  you  would  first  calculate 
their  GSH  level,  use  it  in  the  above  equation  to  get  the  Probability  of  Fatality,  and  then  compare 
that  value  to  a  random  draw.  For  example,  if  a  person  were  to  calculate  a  Probability  of  Fatality 
of  0.5,  and  random  draw  greater  than  0.5,  would  result  in  that  person  becoming  a  fatality. 
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Section  7. 


Vesicant  Countermeasure  Model  Framework 

We  are  including  the  same  medical  countermeasures  in  the  GVM  that  were  proposed  by 
USAMRICD,  with  the  exception  of  calcium  modulators.  As  mentioned  earlier  in  this  paper,  we 
are  excluding  calcium  modulators  because  they  have  not  been  shown  to  have  a  therapeutic  effect 
in  vitro.  The  countenneasures  proposed  by  USAMRICD  are  listed  in  Table  7-1. 


Table  7-1.  Proposed  Medical  Countermeasure  Development  Areas  (Smith  W.  J.,  2009) 


Biochemical  Event 

Pharmacologic  Strategy 

Example 

DNA  alkylation 

Scavengers 

N-acetyl  cysteine 

DNA  strand  breaks 

Cell  cycle  inhibitors 

Mimosine 

PARP  activation 

PARP  inhibitors 

Niacinamide 

Proteolytic  activation 

Protease  inhibitors 

AEBSF  is  a  sulfonyl  fluoride 
compound 

Inflammation 

Anti-inflammatories 

Indomethacin,  capsaicin, 
hydrocortisone 

The  five  countermeasures  we  will  focus  on  have  the  following  mechanisms  of  action: 

1 .  Scavengers  -  Act  as  a  sink  for  vesicants  before  DNA  strand  breakage  can  occur  (Bobb, 
2005).  Used  as  a  prophylactic,  an  intracellular  scavenger  such  as  NAC  increases  the  level  of 
gluthathione  (GSH)  within  cells,  a  compound  that  HD  has  a  natural  affinity  for.  HD  reacts 
with  the  GSH  and  is  rendered  inert.  This  effect  appeared  in  cancer  patients  that  use  HD  as 
chemotherapy.  When  used  in  conjunction  with  NAC,  low  levels  of  HD  have  a  diminished 
effect.  This  prompted  research  into  NAC  as  a  chemoprotectant  against  HD  and  it  has  been 
found  to  significantly  reduce  lung  injury  in  rats,  when  used  as  both  a  prophylactic  and  post¬ 
exposure  countenneasure. 

2.  Cell  cycle  inhibitors  -  Stops  the  proliferation  of  cells  by  stopping  DNA  replication  (Krude, 
1999).  Mimosine,  for  example,  stops  the  replication  of  DNA  in  the  G1  phase,  before  the 
strands  fork.  This  treatment  delays  the  production  of  new  cells,  allowing  DNA  repair  to  take 
place,  and  also  stops  injured  cells  from  attempting  to  replicate,  resulting  in  a  killed  cell. 

3.  PARP  inhibitors  -  Inhibits  poly(ADP-ribose)  polymerase  (PARP),  an  enzyme  responsible 
for  DNA  repair.  It  was  observed  that  PARP  inhibitors  are  capable  of  arresting  DNA 
fragmentation  that  is  associated  with  mustard  injuries  (Meier,  1998). 

4.  Protease  inhibitors  -  Blocks  protease  from  breaking  down  the  epithelial  attachment. 
AEBSF  specifically  targets  proteases  that  cut  peptide  bonds  in  proteins.  The  inhibition 
process  using  AEBSF  is  irreversible,  allowing  DNA  to  repair  itself  and  return  the  body  to  a 
stable  baseline. 
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5.  Anti-inflammatories  -  Multiple  anti-inflammatories  have  been  shown  to  reduce  the  signs 
and  symptoms  associated  with  mustard  exposure.  They  work  to  slow  erythema,  but  their 
effect  decreases  over  time  (Yourick,  1995).  However,  delaying  the  onset  of  erythema  allows 
for  treatment  that  repairs  the  fragmented  DNA. 

In  addition  to  these  five  countermeasures,  we  are  also  including  a  model  for  British  anti- 
Lewisite  (BAL).  BAL  is  a  chelating  agent  that  is  very  effective  at  removing  lewisite  (and  other 
heavy  metals)  at  the  injury  site  (Vilensky  &  Redman,  2003)  and  systemically. 


Figure  7-1.  System  Diagram  for  Vesicant  Countermeasures 


This  section  will  document  the  models  and  data  for  each  of  the  countermeasures.  An 
additional  benefit  to  the  revised  injury  model  mechanism  for  the  GVM  is  we  no  longer  have  to 
derive  rates  of  action  for  each  of  these  countermeasures.  In  all  cases,  the  countermeasures  work 
to  reduce  the  concentration  of  the  compartment  that  it  is  acting  on.  Since  these  individual 
compartments  are  no  longer  based  on  differential  equations,  the  countermeasures  can  now  act  as 
dose-dependent  scaling  factors  that  work  to  reduce  the  concentration  in  the  compartment.  One 
exception  to  this  rule  is  the  cell-cycle  inhibitors,  which  will  effectively  reduce  the  level  of 
mitotic  cell  division  to  zero.  Each  of  the  proceeding  sections  will  discuss  how  the 
countenneasure  model  is  different  than  the  model  proposed  for  HD  and  provide  new  model 
equations  based  on  the  reviewed  literature. 


58 


7.1  Scavengers 


In  the  HD  model,  we  derived  two  rates  that  described  the  efficacy  of  scavengers  against  HD 
exposure:  the  rate  at  which  scavengers  sequester  HD  in  the  lungs,  the  rate  at  which  scavengers 
sequester  HD  in  the  skin.  We  based  our  derivations  on  a  study  performed  by  K.  Atkins,  et.  al. 
(Atkins,  Lodhim,  Hurley,  &  Hinshaw,  2000)  that  studied  the  efficacy  of  N- acetylcysteine  (NAC) 
against  HD  exposures  in  epithelial  cells.  We  will  turn  to  this  study  again  for  discussion  on  how 
to  incorporate  scavengers  into  the  model.  Table  7-2  summarizes  the  result  from  the  Atkins  et  al. 
study. 


Table  7-2.  Efficacy  of  IntracellularSscavenger  NAC  against  HD  (Atkins,  2000) 


Control 

500  pM  HD 

% 

Viable 

% 

Apoptotic 

% 

Necrotic 

% 

Viable 

% 

Apoptotic 

% 

Necrotic 

Untreated 

90.3±3.1 

3.6±2.1 

6.1±2.8 

50.4±4.0 

33.1±5.0 

16.3±4.5 

0.2  inM  BSO 

91.0±3.0 

2.7±0.6 

6.3±3.1 

35.7±8.4 

32.7±10.3 

31.7±17.8 

50  m\I  NAC 

92.3±3.6 

2.5±1.0 

4.3±2.6 

69.3±1.7 

21.2±4.3 

9.5±4.4 

BSO  +  NAC 

92.7±4.0 

2.0±1.0 

5.3±3.2 

51.7±4.7 

37.0±5.7 

11.3±1.2 

2.5  inM  GSH 

88±0.0 

3.6±2.1 

6.1±2.8 

72.7±5.9 

15.7±5.5 

1 1.7±4.2 

5.0  inM  GSH 

88.0±5.5 

4.8±2.5 

9.8±5.6 

75.4±5.6 

11.3±1.3 

1 1.3±4.2 

10.0  inM 
GSH 

83.0±4.0 

4.0±1.7 

13.3±3.0 

66.7±6.4 

10.7±3.2 

23.0±3.6 

The  study  measured  the  viability  of  exposed  cells  (a  value  that  correlates  directly  with  the 
percentage  of  nonnal,  healthy  cells),  and  the  percentages  of  apoptotic  and  necrotic  cells 
following  treatment.  For  a  500pM  exposure  of  HD,  untreated  cells  had  50.4%  viability  with 
33.1%  of  cells  being  apoptotic  and  16.3%  of  cells  being  necrotic.  50mM  of  NAC  increased  the 
percentage  of  viable  cells  to  69.3%  while  decreasing  the  percentage  of  apoptotic  and  necrotic 
cells  to  21.2%  and  9.5%,  respectively.  This  is  equivalent  to  a  protection  factor  of  approximately 
1.5  for  a  1:10  dosage  of  NAC . 

However,  a  second  study  (Sharma,  Vijayaraghavan,  &  Agrawal,  2010)  did  not  result  in  the 
same  level  of  efficacy  for  NAC.  This  study  looked  at  NAC  efficacy  (along  with  other  GSH- 
producing  countenneasures)  in  mice  exposed  to  a  LD50  of  nitrogen  mustard  and  sulfur  mustard. 
In  that  study,  a  dose  of  185  mg/kg  of  NAC  was  applied  following  an  8.1  mg/kg  dosing  of  HD  (at 
comparable  dosage  of  nitrogen  mustard)  and  showed  no  statistically  significant  change  in  the 
percentage  of  viable  cells  at  the  injury  site. 

Obviously,  the  results  from  these  two  studies  are  drastically  different.  One  potential  reason 
for  the  discrepancy  is  the  efficacious  benefit  of  NAC  was  overpowered  by  the  high  dosing  of 
vesicant.  If  we  consider  that  NAC  stimulates  GSH  production  in  the  cell,  it  is  possible  to 
postulate  a  point  at  which  even  the  maximum  rate  of  increased  GSH  production  is  insignificant, 
compared  to  the  dose  of  vesicant  at  the  cellular  level. 
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A  second  possibility  is  that  the  mechanism  of  NAC  administration  in  the  Shanna  et  al.  study 
(oral  dosing)  does  little  to  affect  the  level  of  injury,  but  instead  affects  the  incidence  of  lethality 
by  decreasing  the  level  of  oxidative  stress.  However,  in  the  Shanna  study,  there  was  very  little 
decrease  in  the  GSH  level  in  the  liver  following  NAC  administration. 

A  third  possibility  is  that  NAC  is  not  efficacious  against  percutaneous  exposures.  The  Atkins 
et  al.  study  was  perfonned  on  lung  epithelial  cells  and  the  Shanna  et.  al.  study  was  performed  via 
percutaneous  administration  of  vesicants.  This  possibility  is  strengthened  by  looking  at  a  study 
performed  by  Om  Kumar  et  al.  (Kumar,  Sugendran,  &  Vijayarghavan,  2001).  In  this  study,  the 
researchers  looked  at  the  effect  that  three  scavengers  (trolox,  quercetin,  and  pure  GSH)  have  on 
lethality  in  HD-exposed  mice.  Unfortunately,  NAC  was  not  included  in  the  study.  All  three 
scavengers  showed  a  dramatic  increase  in  survival  time  for  mice  that  were  exposed  via  inhaled 
HD  (100%  of  mice  died  on  day  7  without  treatment;  12.5%,  12.5%,  and  25%  died  on  day  seven 
following  trolox,  quercetin,  and  GSH,  respectively)  and  no  significant  change  in  survival  time 
when  mice  were  exposed  percutaneously. 

Given  the  result  of  these  three  studies,  we  assume  that  NAC  (and  other  scavengers)  are  most 
efficacious  in  the  lungs.  Knowing  that  scavengers  work  to  reduce  the  dose  of  vesicant  available 
to  attack  the  DNA  of  normal,  healthy  cells,  we  use  NAC  to  modify  the  inhaled  dosage  Nn,i-  In 
the  Atkins  et  al.  study,  a  1:10  dosage,  of  NAC  vs.  agent,  results  in  a  39%  increase  in  viable  cells. 
That  seems  like  a  drastic  increase  in  cell  viability,  but  the  Kumar  et  al.  study  showed  the  level  of 
efficacy  that  scavengers  like  NAC  are  capable  of  achieving.  We  incorporate  NAC  into  the  model 
by  accounting  for  its  dose-reducing  properties  in  the  equation  for  D,. 

d  rTmask  (  7-1 ) 

~  J  C  dt  *  b  —  NACreciUction  *  kfoind  *  *  D; 

Since  we  don’t  have  the  necessary  data  to  support  a  dose-dependent  model  of  NAC  efficacy, 
or  found  a  study  that  sampled  multiple  dosings  of  NAC,  or  a  study  of  NAC  being  used  for 
anything  outside  of  prophylaxis,  we  had  to  assume  that  pre-exposure  administration  of  NAC 
leads  to  a  39%  increase  in  viable  cells  in  the  lungs,  as  shown  in  Atkins,  et.  al.  Therefore,  the 
value  for  NACreduction  is  1.0  if  no  NAC  is  administered  and  1.39  if  NAC  is  administered.  Note 
that  the  proposed  application  of  NAC  against  vesicant  injuries  is  as  a  pre-exposure  prophylaxis, 
so  this  model  is  consistent  with  current  doctrine.  If  NAC  is  ever  proposed  as  a  post-exposure 
antidote  for  vesicant  exposures,  we  will  revisit  this  model.  However,  given  the  rapid  action  of 
vesicants  in  the  lungs  and  the  fact  that  NAC  has  not  been  shown  to  be  beneficial  against 
percutaneous  injuries,  we  find  the  likelihood  if  it  being  designated  an  antidote  relatively  small. 

In  order  to  see  what  NAC  does  in  the  GVM,  we  varied  the  inhaled  dose  of  HD  from  250  mg- 
min/m3  (the  lower  bound  of  the  first  dose  band  that  causes  acute  lethality)  to  250  mg-min/m3 
and  calculated  the  value  for  Nk,i.  The  results  are  shown  in  Table  7-3. 
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Table  7-3.  Fraction  of  Killed  Cells  With  and  Without  NAC 


Dose  (mg/man) 

250 

260 

270 

280 

290 

300 

310 

325 

340 

350 

Nk,i  (without 
NAC) 

0.8 

0.83 

0.87 

0.91 

0.96 

1 

1 

1 

1 

1 

Nk,i  (with  NAC) 

0.56 

0.58 

0.61 

0.63 

0.66 

0.69 

0.72 

0.75 

0.79 

0.83 

For  the  high,  lethality-causing  doses,  NAC  significantly  lowers  the  proportion  of  cells.  As  a 
reminder,  the  value  for  inhalation  lethality  is  anything  over  0.8.  The  Nk,i  range  for  the  most 
severe,  non-lethal  AMedP-8  (C)  dose  band  is  0.45  -  0.8.  With  NAC,  the  model  predicts  that  an 

•5 

individual  would  survive  a  dose  of  up  to  340  mg-min/m  ,  but  they  will  still  exhibit  severe  effects. 

7.2  Cell-Cycle  Inhibitors 

When  we  developed  the  cell-cycle  inhibitor  model  for  HD,  we  underwent  a  thorough 
literature  review  and  developed  a  pharmacokinetic  and  dynamic  model  for  cell-cycle  inhibitors. 
The  models  ended  up  predicting  that  cell-cycle  inhibitors  did  more  harm  than  good  as  the  injury 
progressed,  despite  published  studies  showing  a  slight  overall  beneficial  effect  of  cell-cycle 
inhibitors.  At  the  time,  we  reported  that  the  discrepancy  could  come  from  a  lack  of  modeling  for 
apoptotic  vs.  necrotic  death.  Since  cell-cycle  inhibitors  work  to  slow  mitotic  division,  they  would 
work  to  stop  the  apoptotic  pathway  in  the  percutaneous  model.  For  low  doses,  this  would  result 
in  a  lower  concentration  of  apoptotic  cells,  but  at  the  risk  of  having  a  higher  proportion  of 
necrotic  cells.  So,  even  if  the  overall  proportion  of  killed  cells  decreases,  cell-cycle  inhibitors 
would  work  to  actually  increase  the  level  of  necrotic  cells,  thereby  increasing  the  concentration 
we  use  as  a  metric  for  injury  severity. 

We  believe  our  model  captured  the  concentration  and  efficacy  of  cell-cycle  inhibitors,  which 
we  reused  in  the  model  for  the  GVM. 

The  method  for  incorporating  cell-cycle  inhibitors  into  the  HD  model  was  to  multiply  the 
terms  that  account  for  mitotic  cell  division,  by  a  scaling  function  f(CC),  which  is  equal  to: 

f(rr ,  _  0  ,CC<CCdose  (7-2) 

^LL)~  cceff,cc  >  ccdose 

Where  CCdose  is  the  efficacious  dose  of  the  cell-cycle  inhibitor  (for  mimosine  we  derived  a 
value  of  0. 1  ng/mL)  and  Ceff  represented  the  efficacy  of  the  cell-cycle  inhibitor  (for  mimosine  we 
derived  a  value  of  1.0,  meaning  it  was  completely  effective  in  stopping  mitotic  division).  CC  was 
the  calculated  concentration  of  cell-cycle  inhibitor  in  ng/mL,  which  is  given  by: 

CC  =  Dcc  *  y  *  e~zt  (  ) 
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Where  DCC  is  the  user-inputted  administered  dose,  and  y  and  z  were  constants  fit  to  20  and 
0.49,  respectively.  Since  cell-cycle  inhibitors  are  to  be  used  sparingly,  due  to  potential  adverse 
side  effects,  we  modified  the  above  equations  to  account  for  time  of  administration  as  well  as 
administered  dose. 


cc  =  Dcc  *  y  *  e  zt,t>  tadmin  (  7 -4) 

0<  t  <  t-admin 

Where  tadmin  is  the  time  the  cell-cycle  inhibitor  was  given  in  hours. 

The  model  accounts  for  the  effect  of  cell-cycle  inhibitors  by  multiplying  the  mitotic  cell 
division  terms  in  the  healthy  cell  and  injured  cell  equations  by  l-f(CC).  The  relative 
concentration  of  normal,  healthy  cells  therefore  becomes: 

—  NN  =  (l—  /(CC))  *  kSpiit  *  (1  —  Nn )  *  Nn  —  kbind  *  D  *  Nn  +  krepair  *  P  ARP  *  NAD  *  Nt 

And  the  relative  concentration  of  injured  cells  becomes: 

— //  =  (l  —  /(CC))  *  kbind  *  D  *  Nn  —  kspiit  *  (1  —  Nj)  *  Nj  —  krepair  *  P ARP  *  NAD  *  IV,  —  kkm 
*  Prot  *  Ca*  Inflam  *  Nj 

Since  CCeff  is  equal  to  1,  the  implementation  effectively  states  that  the  cell-cycle  inhibitors 
are  completely  effective  as  long  as  they  are  present  in  the  body  above  some  threshold. 

The  model  now  more  accurately  predicts  the  therapeutic  effects  of  cell-cycle  inhibitors 
against  vesicants.  Figure  7-2  shows  the  model  prediction  for  cell-cycle  inhibitor  effectiveness 
against  the  concentrations  of  necrotic  (blue)  and  apoptotic  (red)  cells  as  well  as  the  overall 
concentration  of  killed  cells  (green)  with  and  without  (dashes)  lOOmg  of  mimosine.  As  the  figure 
shows,  cell-cycle  inhibitors  have  a  very  small  effect  on  the  overall  level  of  killed  cells, 
decreasing  the  proportion  of  apoptotic  cells,  and  increasing  the  proportion  of  necrotic  cells. 

The  model  even  does  a  good  job  of  matching  the  published  protection  factors  of  mimosine. 

In  a  study  perfonned  by  ICD  (Smith  W.  S.,  1993),  the  researchers  reported  a  4%  increase  in 
viable  cells  at  24  hours  following  pre-exposure  mimosine  dosing,  whereas  our  model  predicts  a 
3.3%  decrease  in  killed  cells  at  24  hours  (0.152  without  mimosine  and  0.147  with  mimosine); 
however,  the  study  does  not  differentiate  between  apoptotic  and  necrotic  deaths. 
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Figure  7-2.  Cell-Cycle  Inhibitor  Effect  on  Nectrotic  and  Apoptotic  Cells  Following 

Percutaneous  Exposure  to  Vesicants 

Figure  7-3  shows  the  model  prediction  for  cell-cycle  inhibitor  effectiveness  for  five  different 
dosing  times.  Notice  how  the  cell-cycle  inhibitor  is  predicted  to  be  more  efficacious  as  the 
dosing  time  gets  later.  This  agrees  with  the  model’s  framework,  as  cell-cycle  inhibitors  provide  a 
window  of  protection  against  apoptotic  cell  death.  As  that  window  is  shifted  towards  the  time  at 
which  injury  levels  are  at  their  highest  (12-48  hours  according  to  the  model),  the  cell-cycle 
inhibitor  would  have  the  greatest  effect  against  cells  dying  via  attempting  mitotic  division.  As 
the  dosing  time  is  delayed  further,  however,  the  cell-cycle  inhibitors  will  begin  to  lose  efficacy 
according  to  our  model.  This  trend  is  also  noted  in  the  literature.  In  the  Smith  et.  al.  study, 
mentioned  above,  the  efficacy  of  cell-cycle  inhibitors  against  HD  was  measured  when  mimosine 
was  dosed  at  3,  5,  and  21 -hours  post-exposure.  The  efficacy  at  24  hours  for  the  three  dosings  was 
measured  to  be  17%,  15%,  and  15%  respectively,  compared  to  the  4%  when  the  cells  were  dosed 
before  exposure. 
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Figure  7-3.  Cell-Cycle  Inhibitor  Effectiveness  for  Different  Administration  Times 
7.3  PARP  Inhibitors 

For  the  FID  model,  PARP  inhibitors  were  measured  by  the  reduction  of  the  concentration  of 
PARP.  We  defined  an  inhibition  rate  kPARPi  that  permanently  removed  PARP  from  the  model. 
Once  again,  the  model  produced  results  that  were  not  supported  by  the  literature.  In  one  study 
(Bhat,  Ramachandra,  J.,  &  Ray,  2004),  the  researchers  found  that  PARP  inhibition  using  3- 
amino  benzamide  (3-AB)  reduced  the  repair  rate  of  human  epidermal  keratinocytes  by  66%. 
Whereas  the  HD  model  did  predict  a  reduction  in  the  concentration  of  killed  cells  early  on  (<15 
hours),  the  PARP  inhibitors  ultimately  resulted  in  an  increase  in  cellular  death  in  our  model 
when  compared  to  untreated  individuals.  Looking  back  at  the  model,  this  probably  resulted  from 
not  having  a  toxicodynamic  model  for  vesicants  that  accounted  for  the  rate  at  which  HD  is 
metabolized.  Nor  did  we  have  a  model  that  differentiated  between  necrotic  and  apoptotic  cell 
death.  With  the  GVM,  we  are  now  accounting  for  PARP  inhibitors  by  first  including  a 
pharmacokinetic  model  for  PARP  and  then  using  the  concentration  of  PARP  inhibitors  to  scale 
the  concentration  of  PARP  in  the  body. 

For  the  pharmacokinetic  model  we  turn  to  a  study  conducted  by  Rolf  W.  Sparidans 
(Sparidans,  Martensa,  Iersel,  Hartigh,  Schellensa,  &  Beijnena,  2011)  where  plasma 
concentrations  of  the  PARP  inhibitor  olaparib  were  measured  in  a  human,  following  a  50mg 
dose.  The  result  for  this  study  is  shown  in  Figure  7-4  along  with  the  corresponding  quadratic  fit. 
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Olaparib  Concentration  in  the  Plasma  as  a 
Function  of  Time 


Figure  7-4.  PARP  Concentrations  in  the  Plasma  Following  a  50mg  dose 

The  PARP  inhibitor  concentrations  remain  relatively  flat  for  the  first  six  hours  before 
experiencing  a  rapid  drop-off  in  concentration.  The  equation  for  PARP  inhibitor  concentration 
(in  ng/mL)  is  given  by: 


PARPIn  =  -4.6  *  (t  -  tadmin)2  +  28.2  *  (t  -  tadmin )  +  87 


(7-5) 


where  tadmin  is  the  time  at  which  the  PARP  inhibitor  is  given.  Using  this  model,  an  individual 
will  get  approximately  8  hours  of  protection  from  a  single  dosing  of  PARP  inhibitors.  We  next 
defined  a  parameter  PARPIneff,  which  represents  the  efficacy  of  PARP  inhibitors.  Using  the 
study  perfonned  by  Bhat  et  ah,  we  set  the  value  of  PARPIneff,  so  that  there  would  be  a  66% 
reduction  in  PARP  levels  when  PARP  inhibitors  are  at  their  highest  concentrations.  We  do  this 
by  first  nonnalizing  Equation  6-5  to  have  a  maximum  value  of  1.0,  as  shown  in  Figure  7-5. 

This  modifies  Equation  (  7-5  )  to: 


PARPIn  =  -0.03  *  (t  -  tadmin)2  +  0-19  *  (t  -  tadmin)  +  0-58 


(7-6) 
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Olaparib  Concentration  (Normalized  to  1)  in 
the  Plasma  as  a  Function  of  Time 
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Figure  7-5.  PARP  Concentrations  (Normalized  to  1)  in  the  Plasma  Following  a  50mg  Dose 


The  concentration  of  PARP  can  now  be  modeled  as: 

PARP  =  (1  -  P ARPeff  *  PARPconc )  *  5.91  *  Nlx  (  7'7> 

The  model  prediction  for  PARP  inhibitor  effectiveness  for  three  different  dosing  times  is 
shown  in  Figure  7-6.  Notice  that  PARP  inhibitors  have  the  opposite  effect  of  cell-cycle 
inhibitors;  they  actually  work  to  reduce  the  percentage  of  necrotic  cell  death,  while  increasing 
the  percentage  of  apoptotic  cells.  This  is  consistent  with  the  model  and  with  the  known  injury 
mechanism  for  vesicants.  Also,  PARP  inhibitors  are  predicted  to  be  more  effective  with  delayed 
treatment  times,  as  the  8-hour  window  of  protection  is  most  beneficial  when  PARP  would  be 
most  active.  One  interesting  fact  about  the  model  is  that  if  you  delay  the  dosing  long  enough, 
PARP  inhibitors  actually  increase  the  percentage  of  necrotic  cells.  This  is  probably  due  to 
NAD+  concentrations  that  are  directly  correlated  with  PARP  concentrations  in  this  model. 
Keeping  PARP  levels  low  at  a  time  when  NAD+  levels  are  low  leads  to  an  increase  in  necrotic 
cell  death.  We  have  no  evidence  to  support  this  trend,  so  for  this  implementation  we  recommend 
keeping  the  PARP  inhibitor  dosing  at  less  than  10  hours. 
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Figure  7-6.  PARP  Inhibitor  Effectiveness  for  Three  Different  Dosing  Times 
7.4  Protease  Inhibitors 

We  did  a  comprehensive  review  of  protease  inhibitors  when  developing  the  HD  model,  but 
once  again,  accounted  for  the  effect  of  protease  inhibitors  by  defining  another  set  of  rate 
equations.  We  designed  these  equations  to  irreversibly  lower  the  concentration  of  protease.  For 
the  GVM,  we  can  model  protease  inhibitors  the  same  way  we  model  PARP  inhibitors,  by  first 
starting  with  a  toxicokinetic  model  for  protease  concentration,  normalizing  it  to  1.0,  and  then 
defining  a  maximum  efficacy  for  the  protease  inhibitor.  Unfortunately,  we  have  not  yet  found  a 
study  measuring  the  pharmacokinetics  of  protease  inhibitors  in  the  literature,  so  we  had  to 
assume  that  the  concentration  of  protease  inhibitors  follows  a  similar  trend  as  the  concentration 
of  PARP  inhibitors. 

This  would  make  the  concentration  of  protease  inhibitors  (in  ng/mL)  equal  to: 

Protconc  =  -0.03  *  (t  -  tadmin )2  +  0.19  *  (t  -  tadmin )  +  0.58  (  7~8 

We  will  continue  our  literature  search  and  attempt  to  find  a  pharmacokinetic  study  of 
protease  inhibitors  in  humans,  in  order  to  update  this  portion  of  the  model. 

To  fit  the  maximum  efficacy  of  protease  inhibitors  we  can  turn  to  several  studies.  One  study 
on  protease  inhibitor  efficacy  against  inhaled  HD  in  rats  showed  a  significant  benefit  when 
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aprotinin  was  used  to  irreversibly  bind  with  protease  (Anderson,  2009).  Aprotinin-treated  rats 
showed  a  50-percentage  point  decrease  in  bronchial  injuries  and  a  45-percentage  point  difference 
in  edema,  effectively  eliminating  severe  respiratory  effects  from  the  animals. 

Another  study  conducted  in  vivo  using  a  mouse  ear  model,  tested  the  protective  efficacies  of 
multiple  protease  inhibitors  (Casillas,  Therapeutic  Approaches  to  Dennatotoxicity  by  Sulfur 
Mustard  I.  Modulation  of  Sulfur  Mustard-induced  Cutaneous  Injury  in  the  Mouse  Ear  Vesicant 
Model,  2000).  Of  the  five  protease  inhibitors  tested,  only  one  (peptide  phosphonate)  showed 
significant  reduction  of  edema  (16%  reduction  from  positive  control),  but  none  had  significant 
reduction  in  subepidennal  blisters,  epidennal  necrosis,  or  contralateral  epidermal  necrosis.  One 
point  of  note  is  that  the  protease  inhibitors  were  only  used  as  pretreatment.  Looking  back  at  our 
model  output  (and  the  measured  protease  levels  it  was  based  on)  it  is  apparent  that  protese  levels 
take  several  hours  to  rise.  Several  hours  to  several  days  are  needed  for  dead  cells  to  accumulate, 
which  would  result  in  blistering.  It  is  possible  that  protease  inhibitors  are  simply  not  beneficial  as 
prophylaxis,  much  like  PARP  inhibitors  and  cell-cycle  inhibitors. 

One  other  study  looked  at  the  efficacy  of  using  the  protease  inhibitor  doxycycline  (DOX) 
against  HD  and  using  nitrogen  mustard  against  human  bronchial  epithelial  cells  in  vitro 
(Rappeneau,  2000).  Two  treatment  strategies  were  considered  in  this  study,  one  with  immediate 
co-treatment  and  one  with  2.5  to  90-minutes  delayed  treatment.  For  the  delayed  treatment  model, 
the  researchers  looked  at  treatment  using  NAC  alone  and  NAC  in  combination  with  DOX.  The 
benefits  of  NAC  alone  against  HD  were  mainly  seen  in  the  first  5-10  minutes,  with  a  peak 
protection  factor  of  1.6  taking  place  2.5  minutes  post-exposure.  Also,  DOX  on  its  own  had  a 
protection  factor  of  1.3.  However,  NAC  combined  with  DOX  offered  a  protection  factor  of  2.0 
that  peaked  at  2.5-minutes  post-exposure,  but  there  was  still  significant  protection  with  the 
combination  treatment  if  administered  within  60-minutes. 

We  decided  to  use  the  Rappeneau,  et.  al.  study,  as  it  was  conducted  on  human  cell  cultures. 
Additionally,  the  toxicokinetic  model  will  handle  the  time-dependent  differences  in  efficacy.  A 
peak  protection  factor  of  1 .3  is  obtained  using  Berkeley  Madonna  to  derive  a  value  for  Proteff 
used  in  the  equation  below: 

Protease  =  (1  —  Proteff  *  Protconc) 35. 6  *  e~3A*ATP  (  7-9  ) 

We  need  to  make  the  comparison  at  the  time  of  peak  protease  inhibitor  concentration  in  order 
to  correctly  use  the  in  vitro  data.  If  protease  inhibitors  are  administered  at  the  time  of  exposure, 
the  peak  concentration  time  is  at  approximately  three  hours.  The  value  for  Pro teff  that  results  in  a 
33%  decrease  in  necrotic  cells  at  three  hours  is  0.33. 
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Figure  7-7  shows  the  concentration  of  necrotic  cells  using  three  different  administration 
times  of  protease  inhibitors.  The  model  predicts  an  initial  benefit  from  protease  inhibitors  but 
little  efficacy  as  time  goes  on.  At  this  point,  we  are  going  to  use  the  model  to  predict  the  efficacy 
of  continuous  dosings  of  protease  inhibitors  every  12  hours,  for  a  48-hour  period. 


Protease  Inhibitor  Effectiveness  for  Different 
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Figure  7-7.  Protease  Inhibitor  Effectiveness  for  Three  Different  Treatment  Times 
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Figure  7-8  shows  the  effectiveness  of  protease  inhibitors  when  given  regularly  over  a  48- 
hour  period.  Although  we  did  see  an  increased  efficacy  prediction  from  the  model,  the  benefit 
was  not  substantial.  However,  we  do  not  have  data  to  verify  whether  or  not  this  is  a  valid 
conclusion. 


Protease  Inhibitor  Effectiveness  When  Dosed 
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Figure  7-8.  Protease  Inhibitor  Effectiveness  for  Multiple  Dosings 
7.5  Anti-inflammatories 

For  the  HD  model,  we  developed  a  pharmacokinetic  and  pharmacodynamic  model  of  anti¬ 
inflammatory  effectiveness  using  the  published  literature.  There  is  a  significant  amount  of 
literature  on  anti-inflammatory  effectiveness  as  it’s  the  only  countermeasure  being  considered 
that  was  actually  used  against  human  vesicant  exposures.  There  are  two  significant  departures 
from  the  HD  model,  however.  First,  there’s  the  fact  that  we  used  the  anti-inflammatory 
toxicokinetic  model  to  modify  the  rate  equations  for  anti-inflammatories.  Once  again,  we  are  no 
longer  using  rate  equations  for  the  underlying  biological  processes.  Secondly,  the  inflammatory 
response  model  developed  for  HD  was  not  based  on  any  published  data,  a  limitation  that  we 
corrected  in  the  GVM.  However,  we  reused  the  anti-inflammatory  model  that  was  developed, 
since  the  efficacy  and  concentration  profiles  for  anti-inflammatories  were  well  founded;  we  did 
change  the  implementation. 
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The  equation  for  anti-inflammatory  efficacy  defines  a  rate  of  change  in  anti-inflammatory 
concentration  (x).  Using  these  parameters  and  an  input  dose  (A Iinlamd0se)  in  mg,  we  defined  the 
equation  for  the  concentration  of  effective  anti-inflammatories  to  be: 

Ain  flam  conc  =  AInflamdose  *  e ~x*1  (  7~10  ) 

For  the  GVM  we  will  modify  that  equation  to  account  for  administration  time  (tadmin),  as 
shown  below: 


AInflamconc  =  AIinflamdose  *  e  x*(f  tadmin)  (7-11) 

In  the  HD  model,  we  also  derived  parameter  values  for  the  anti-inflammatories  olvanil  and 
indomethacin,  as  shown  in  Table  7-4. 

Table  7-4.  Rate  of  Change  for  the  Anti-Inflammatories  Olvanil  and  Indomethacin 


X 

Olvanil 

.1 

Indomethacin 

.1 

In  the  HD  model  we  then  used  these  values  to  define  a  rate  of  anti-inflammatory 
effectiveness,  kinflam-  However,  we  will  use  the  anti-inflammatory  concentration  in  a  manner 
similar  to  PARP  and  protease  inhibitors  by  first  nonnalizing  the  concentration  to  1  and  getting 
the  following  equation: 


AInflamconc  =  e  x*<~t  1  admin )  (7-12) 

We  then  modified  the  equation  for  inflammatory  cell  concentration  to  get: 

Inflam  =  (l  —  AInflameff  *  Inflamconc )  *  0.98  *  exp(Nf )  C  7-H  ) 

where  AInflamefi  is  the  efficacy  of  the  anti-inflammatory.  We  then  revisit  the  study 
perfonned  by  Robert  Casillas  et  al.  (Casillas,  2000)  that  was  used  to  define  the  rate  parameters  in 
the  HD  model,  to  define  the  value  of  PJnflameff.  Casillas  showed  that  there  was  a  significant 
reduction  in  edema  (19%  and  23%  at  24  hours  post-exposure  for  olvanil  and  indomethacin, 
respectively)  when  the  anti-inflammatories  were  applied  topically  at  10  minutes  post-exposure. 
Since  edema  is  a  function  of  necrotic  cell  death  in  our  percutaneous  model,  we  can  use  these 
results  to  fit  a  value  for  Ieff  in  the  GVM.  The  resulting  values  are  shown  in  Table  7-5. 


Table  7-5.  Efficacies  for  the  Anti-Inflammatories  Olvanil  and  Indomethacin 


AInflameff 

Olvanil 

1.1 

Indomethacin 

1.3 
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The  model’s  prediction  of  anti-inflammatory  effectiveness,  administering  a  dosing  every  12 
hours  for  48  hours  is  shown  in  Figure  7-9.  The  model  predicts  a  significant  decrease  in  necrotic 
cell  death,  following  anti-inflammatory  application.  This  decrease  is  larger  than  we’ve  seen  with 
other  countermeasures  included  in  this  model,  except  for  scavengers.  Continued  application  of 
anti-inflammatories  would  lead  to  an  even  greater  reduction  in  necrotic  cells,  unlike  other 
treatments  included  in  this  model. 
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Figure  7-9.  Anti-Inflammatory  Effectiveness  for  Multiple  Dosings 
7.6  British  Anti-Lewisite  (BAL) 

BAL  is  an  incredibly  efficacious  heavy-metal  chelating  agent  that  works  very  well  against 
arsenic  compounds,  like  lewisite.  We  have  not  found  any  modern  studies  that  measure  the 
effectiveness  of  BAL  against  lewisite,  but  that  is  likely  because  early  studies  have  proven  that 
the  countermeasure  is  nearly  100%  effective  if  administered  early.  One  experiment  in  rabbits 
showed  a  near  complete  recovery  using  a  5%  solution  of  BAL,  against  enough  lewisite  to  destroy 
the  eyes  of  rabbits.  Similar  experiments  were  conducted  on  rabbits  with  percutaneous  exposures 
to  lewisite  and  even  on  human  volunteers  (Waters  &  Stock,  1945).  The  only  downside  to  BAL  is 
that  nearly  50%  of  patients  exhibit  side  effects  after  receiving  a  5mg/kg  dose  intramuscularly 
(Klassen,  2002).  In  34  humans  given  intramuscular  BAL  doses  of  3.6  to  5.5  mg/kg,  there  were 
some  observed  side  effects  (mainly  headaches  and  nausea),  but  all  subjects  recovered  after  two 
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hours  with  no  cumulative  toxicity  by  multiple  injections  (Sulzberger,  Baer,  &  Kanof,  1946).  We 
have  not  found  any  evidence  to  suggest  that  BAL  is  not  100%  effective  at  reducing  the  dose  of 
lewisite  from  the  time  of  administration.  Nor  will  continued  administration  will  lead  to  anything 
outside  of  mild  side  effects,  most  of  which  are  likely  not  even  operationally  relevant. 

Due  to  these  assumptions,  our  model  for  BAL  is  included  in  the  following  equations: 


Do  —  0,t  >  BALadmin 

(7-14) 

Di  —  0,t  >  BALadmin 

(7-15) 

Dp  —  0 ,t>  BALadmin 

(7-16) 

Ds  —  0,  t  >  BALadmin 

(7-17) 

To  conclude,  BAL  completely  removes  the  dose  from  the  eyes,  lungs,  skin,  and  systemic 
system  from  the  moment  it  is  administered. 
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Section  8. 


Example  Problem 

For  illustrative  purposes  we  now  present  a  walk-through  of  the  model.  The  scenario  is  as 
follows: 

An  individual  is  exposed  to  FID  via  a  line-source  located  upwind.  The  attack  sets  off  a 
detector  and  the  alarm  goes  out  within  seconds.  The  individual,  however,  was  instantly  exposed 
to  100  mg  of  liquid  HD  and  a  10  mg/m3  concentration  of  HD  vapor  for  five  minutes,  before 
putting  on  his  mask.  Two  minutes  later,  he  is  in  full  MOPP  gear.  The  individual  does  not 
decontaminate  his  skin  or  his  eyes.  The  resulting  injury  to  the  eyes,  lungs,  and  skin  are  shown  in 
Figure  8-1. 


Figure  8-1.  Eye,  Lung,  and  Skin  Injury  Levels  and  Killed  Cell  Concentrations  for 

Combined  Exposure 


The  model  predicts  that  the  individual  will  be  in  the  middle  of  the  lowest  percutaneous  dose 
band  from  AMedP-8(C),  on  the  border  of  the  2nd  and  3ld  ocular  dose  band,  and  on  the  border  of 
the  lowest  inhalation  dose  band.  This  means  that  the  individual  will  reach  Severity  Level  1,  for 
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skin  effects  in  36  hours;  reach  Severity  Level  1  in  four  hours  and  Severity  Level  2  in  six  hours, 
for  ocular  effects,  and  no  respiratory  symptoms;  but,  the  individual  would  reach  Severity  Level  1 
for  upper  GI  effects  at  8  hours.  The  GSH  suppression  level  is  so  low  that  there  is  no  risk  of 
becoming  a  fatality. 

We  assume  that  the  individual  did  not  have  NAC  before  exposure,  so  there  is  no  therapeutic 
effect  from  bioscavengers.  However,  if  the  individual  was  given  NAC  before  exposure,  the 
concentration  of  killed  cells  in  the  lungs  would  have  been  reduced  from  0.16  to  0.12,  no  longer 
making  him  an  inhalation  casualty. 

At  one  hour,  the  individual  receives:  100  mg  of  cell-cycle  inhibitors,  a  single  dosing  of 
PARP  inhibitors,  protease  inhibitors,  and  anti-inflammatories.  The  cell-cycle  inhibitors  have  a 
very  minor  effect  on  the  concentration  of  killed  cells  in  the  eye  (the  value  drops  from  0.38  to 
0.36),  but  this  is  enough  to  put  him  in  the  lowest  dosage  band  for  ocular  injuries,  avoiding  a 
Severity  Level  2  injury.  The  combination  of  medications  have  a  greater  effect  on  the  skin, 
reducing  the  concentration  of  necrotic  cells  from  0.3 1  to  0.23,  but  this  decrease  is  not  enough  to 
move  the  individual  out  of  the  first  dose  band  and  we  will  still  have  mild  skin  effects. 

The  same  type  of  analysis  can  be  done  for  different  exposures  and  combinations  of 
medications,  but  the  end  result  should  always  relate  to  the  AMedP-8(C)  dosage  bands  and  the 
lethality  equation  for  systemic  dose. 
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Section  9. 


Summary  of  Assumptions  and  Limitations 

The  following  is  a  summary  of  assumptions  and  limitations  of  the  GVM  and  its  underlying 

data  sources. 

1 .  Exposure 

2.  The  model  uses  the  same  percutaneous  toxicokinetic  model  for  all  vesicants.  However, 
different  vesicants  will  have  different  rates  of  evaporation  and  diffusion  that  should  be 
addressed  by  the  model.  We  will  aim  to  account  for  these  differences  in  the  third  year,  as  the 
models  undergo  verification  and  validation  and  are  transitioned  into  useable  source  code. 

3.  We  have  not  yet  implemented  a  toxicokinetic  model  for  inhalation  and  ocular  exposures. 

4.  The  model  does  not  consider  ingestion  or  intravenous  exposure. 

5.  The  only  sink  for  vesicants  in  the  lungs  and  eyes  is  healthy  cells.  We  do  not  consider  GSH 
suppression  or  other  potential  metabolic  routes  for  these  two  models. 

6.  We  assume  that  the  three  nitrogen  mustards  have  the  same  toxicity.  We  assume  that  lewisite 
has  the  same  toxicity  as  HD. 

7.  Injury  and  lethality 

8.  We  assume  that  PARP  concentration  is  solely  dependent  on  the  concentration  of  injured 
cells. 

9.  We  assume  that  NAD+  concentration  is  solely  dependent  on  the  concentration  of  PARP. 

10.  We  assume  that  that  ATP  concentration  is  solely  dependent  on  the  concentration  of  NAD+. 

1 1.  We  assume  that  protease  concentration  is  solely  dependent  on  the  concentration  of  ATP. 

12.  We  assume  that  inflammatory  response  is  solely  dependent  on  the  concentration  of  injured 
cells. 

13.  We  assume  that  apoptotic  cell  death  is  solely  dependent  on  mitotic  division.  There  is 
evidence  that  this  may  not  be  the  case,  but  we  have  yet  to  find  data  on  which  to  base  a 
different  mechanism. 

14.  We  assume  that  necrotic  cell  death  is  solely  dependent  on  the  depletion  of  ATP. 

15.  We  assume  that  the  concentration  of  necrotic  cells  is  an  indicator  of  severity  of  injury. 

16.  We  assume  that  necrotic  cell  concentration  can  be  used  to  detennine  inhalation  fatalities. 

17.  We  assume  that  GSH  suppression  can  be  used  to  determine  systemic  fatalities. 

18.  Treatment 

19.  We  have  not  yet  developed  a  true  pharmacokinetic  model  for  protease  inhibitors.  We  used 
the  PARP  inhibitor  model.  This  is  a  limitation  that  we  would  like  to  revisit,  as  we  begin 
implementation  of  the  models. 
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20.  We  assume  that  BAL  is  100%  effective,  from  the  moment  it  is  administered.  A 
phannacokinetic  model  of  BAL  would  likely  show  a  difference  in  BAL  effectiveness  as  a 
function  of  time. 

2 1 .  Several  data  sources  for  countermeasure  effectiveness  come  from  in  vivo  experiments.  The 
model  should  be  calibrated  to  in  vitro  data,  when  available.  We  will  continue  our  literature 
review  to  see  if  there  is  better  data  available  for  the  models. 
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Section  10. 


Conclusions  and  Recommendations 

When  we  set  out  to  create  a  generic  vesicant  model,  we  assumed  that  we  could  start  from  the 
HD  framework  and  refit  all  of  the  parameters  for  each  vesicant  of  interest.  We  quickly  realized 
that  we  would  be  doubling  our  parameter  space  with  each  new  agent  and  that  we  would  not  have 
a  truly  “generic”  model  of  vesicant  exposures.  In  order  to  truly  make  the  model  generic,  we 
needed  to  revisit  all  of  the  literature  and  decouple  the  biological  processes  from  the  effecting 
agent.  This  paper  represents  all  of  the  work  that  went  into  that  process. 

We  believe  we  have  a  better  model,  as  a  result  of  the  redesign.  All  of  the  biological  process 
models  are  based  on  published  studies  using  human  cell  exposures.  We  reduced  our  rate 
equations  down  to  four  and  our  rate  parameters  down  to  three,  with  only  one  of  them  being 
intrinsically  tied  to  the  vesicant. 

While  improving  on  the  model’s  design,  we  added  several  sub-models  to  the  injury 
mechanism  that  allowed  us  to  better  quantify  the  effect  of  medical  countermeasures.  We  added  a 
toxicokinetic  model  for  percutaneous  exposure  that  accounts  for  vesicant  diffusion  through  the 
skin,  evaporation  into  the  atmosphere,  and  the  effect  of  metabolizing  agents.  We  included  a 
distinction  between  necrotic  and  apoptotic  cell  death,  to  better  understand  the  effects  of 
countenneasures.  We  included  glutathione  suppression  and  systemic  distribution,  to  account  for 
percutaneous  fatalities. 

We  also  documented  our  equation  development  and  fitting  processes,  so  we  can  quickly 
incorporate  new  data  as  it  becomes  available.  This  model  requires  significant  amounts  of  data. 
We  hope  that  a  continued  interest  in  toxicodynamic  models  will  lead  to  new  studies  aimed  at 
capturing  the  biological  processes,  resulting  from  vesicant  exposure  and  reporting  the  data  in  a 
way  that  is  beneficial  to  the  model. 

The  true  benefit  of  this  model  is  that  it  provides  a  mechanism  for  postulating  the  therapeutic 
effects  of  current,  proposed,  and  future  countenneasures,  aimed  at  combating  vesicant  exposures. 
It  can  and  should  be  refined  in  the  future,  as  we  begin  the  process  of  publishing  the  model  in 
peer-reviewed  literature  and  presenting  it  at  conferences.  We  believe  this  model  represents  a 
great  starting  point  in  modeling  the  effect  of  vesicant  exposures  and  countermeasures;  yet,  there 
remains  room  for  further  development.  There  are  mechanisms  of  injury  that  likely  have  not  yet 
been  discovered,  countermeasures  that  have  not  yet  been  proposed,  and  data  sources  that  have 
not  yet  been  studied.  But,  the  framework  was  created  to  be  scalable,  and  it  is  our  intention  to 
improve  this  model  as  researchers  continue  to  study  the  effects  of  vesicants  and  their 
countenneasures. 
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Section  12. 


Definitions,  Acronyms,  and  Abbreviations 


lxLD50 

A  dose  that  is  1.0  times  the  LD50 

ADP 

Adenosine  Diphosphate 

AEBSF 

4-(2-Aminoethyl)  benzenesulfonyl  fluoride  hydrochloride 

AMedP-8 

Allied  Medical  Publication  8 

ARA 

Applied  Research  Associates,  Inc. 

ATP 

Adenosine  Triphosphate 

ARDS 

Acute  Respiratory  Distress  Syndrome 

BAL 

British  anti-Lewisite 

BAPTA 

1 ,2-bis(o-aminophenoxy)ethane-N,N,N’,N'-tetraacetic  acid 

CB 

Chemical  and  Biological 

CBD 

Chemical  and  Biological  Defense 

CDC 

Centers  for  Disease  Control 

DoD 

Department  of  Defense 

DTRA 

Defense  Threat  Reduction  Agency 

EPA 

Environmental  Protection  Agency 

Gryphon 

Gryphon  Scientific,  LLC 

GSH 

Glutathione 

GVM 

Generic  Vesicant  Model 

HD 

Sulfur  Mustard 

HN-1,  HN-2.  HN-3 

Nitrogen  Mustard  Analogues 

JPEO-CBD 

Joint  Program  Executive  Office  for  Chemical/Biological 
Defense 

JPM-IS 

Joint  Program  Manager-Infonnation  Systems 

JSTO 

Joint  Science  and  Technology  Office 

LD50 

Median  lethal  dose 

LDX 

Dose  causing  X%  mortality 

MCM 

Medical  Countenneasure 

NAG 

N-acetylcysteine 

OP 

Organophosphate 
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NAD+ 

PARP 

PBPK/PD 

RD 

RD-CBI 

SM 

TBD 

TK/TD 

USAMRICD 

USAMRIID 


Nicotinamide  Adenine  Dinucleotide 
Poly(ADP-ribose)  polymerase 

Physiology  Based  Phannacokinetic/Pharmacodynamic 
Research  &  Development 

Infonnation  Systems  Capability  Development  Division 
Sulfur  Mustard 
To  Be  Determined 
T  oxicokinetic/T  oxicodynamic 

US  Army  Medical  Research  Institute  for  Chemical  Defense 
US  Army  Medical  Research  Institute  for  Infectious  Disease 
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Section  13. 


Percutaneous  Source  Code 


METHOD  RK4 

STARTTIME  =  0 
STOPTIME=1400 
DT  =  0.02 

{Exposure  calculations} 

Vapor  =  0  ;Average  concentration  in  mg/m3 

Liquid  =  100  ;Average  concentration  in  mg/m2 

Tmask  =  1/60 

Tmopp=  1/60 

CF  =  1.2 

Sexposed=.93 

BreathingRate=0.9 


Skin=BasalEpidermis+Dermis 

d/dt(Evaporated)  =  k48*Surface 
init  Evaporated  =  0 

d/dt(Surface)  =  IF  TIME  <=  Tmask  THEN  Vapor*CF  -k48*Surface-  k49*Surface+  k94*StratumCorneum 
ELSE  IF  TIME  <=Tmopp  THEN  Vapor*CF*Sexposed  -k48*Surface-  k49*Surface+  k94*StratumCorneum 
ELSE  -k48*Surface-  k49*Surface+  k94*StratumCorneum 
init  Surface  =  Liquid 

d/dt(StratumCorneum)  =  k49*Surface+k29*BasalEpidermis-k92*StratumCorneum-k94*StratumCorneum 
init  StratumCorneum=  0 

d/dt(BasalEpidermis)  =  k92*StratumCorneum+k32*Dermis-k23*BasalEpidermis- 

k29*BasalEpidermis+k22_2*Metabolites-k2_22*BasalEpidermis 

init  BasalEpidermis=  0 

d/dt(Dermis)  =  k23*BasalEpidermis+k1 3*Vascular+k53*Fat-a_k31  *exp(-b_k31  *TIME)*Dermis- 

k32*Dermis-k35*Dermis 

init  Dermis=  0 

d/dt(Fat)  =  k35*Dermis-k53*Fat 
init  Fat=  0 

d/dt(Metabolites)  =  -k22_2*Metabolites+k2_22*BasalEpidermis 
init  Metabolites  =  0 

d/dt(Vascular)=a_k31*exp(-b_k31*TIME)*Dermis-k13*Vascular 
init  Vascuiar=0 

k48=0. 34*60 
k49=0.25*60 
k94=0.23*60 
k92=0.26*60 


89 


k29=0.22*60 
k23=0.24*60 
k32=0. 12*60 
a_k31  =  0.33*60 
b_k31  =0.0001  *60 
k13=0. 12*60 
k35=0. 0058*60 
k53=0. 0060*60 
k2_22  =  0.12*60 
k22  2  =  0.007*60 


d/dt(Systemic)  =  Dermis*(exp(-TIME/taob)-exp(-TIME/taod)) 
init  Systemic  =  0 
taob  =  39 
taod  =  0.7 

;Equations  for  Skin  Injury 

d/dt(NN)  =  k_split*(1-NN)*NN*fcc  -  k_bind*Skin*NN+k_repair*PARP*NAD*NI 
init  NN=1 .0 

Skinlnjury  =  1-NN 

d/dt(NI)  =  k_bind*Skin*NN  -  k_repair*PARP*NAD*NI  -  k_kill*Protease*lnflam*Ca*NI  -  k_split*(1-NI)*NI*fcc 
init  Nl  =  0.0 

CC  =  IF  TIME<  AdminTime  THEN  0  ELSE  Dcc*a*exp(-b*(TIME-AdminTime)) 

Dcc=0 
a  =  20 
b  =  0.49 

fCC  =  IF  CC  <  CCDose  THEN  1  ELSE  0 

CCDose  =  0.1 

AdminTime=1 

d/dt(NKN)  =  k_kill*Protease*lnflam*Ca*NI 
init  NKN  =  0.0 

d/dt(NKA)  =  k_split*(1-NI)*NI*fcc 
init  NKA  =  0.0 

Totallnjury  =  NKA+NKN 

A  =1 .5*NAD 

PARP  =  (1-PARPi*PARPic)*5.9088*NI 

;PARPi  =  .66 

PARPi=0 

PARPic  =  -0.03*(TIME-tparpadmin)A2+.19*(TIME-tparpadmin)+0.58 
limit  PARPic  >=  0 

tparpadmin=0 

NAD  =1.0827*EXP(-.511*PARP) 

Protease  =  (1-Proteff*Protconc)*35.6*exp(-3.4*A) 

Protconcl  =  -0.03*(TIME-tprotadmin1  )A2+.1 9*(TIME-tprotadmin1  )+0.58 
Protconc2=  -0.03*(TIME-tprotadmin2)A2+.19*(TIME-tprotadmin2)+0.58 
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Protconc3  =  -0.03*(TIME-tprotadmin3)A2+.1 9*(TIME-tprotadmin3)+0.58 
Protconc4=  -0.03*(TIME-tprotadmin4)A2+.19*(TIME-tprotadmin4)+0.58 
Protconc5=  -0.03*(TIME-tprotadmin5)A2+.19*(TIME-tprotadmin5)+0.58 
limit  Protconcl  >=  0 
limit  Protconc2  >=  0 

limit  Protconc3  >=  0 
limit  Protconc4  >=  0 
limit  Protconc5  >=  0 

Protconc  =  Protconcl +Protconc2+Protconc3+Protconc4+Protconc5 

tprotadmin1=0 

tprotadmin2=12 

tprotadmin3=24 

tprotadmin4=36 

tprotadmin5=48 

;Proteff=0.33 
Proteff  =  0 

Inflam  =  (1-leff*lconc)*.9811*exp(0.0638*NI) 
limit  Inflam  >=0 

Iconcl  =  exp(-xinflam*TIME-tinflamadmin1) 

lconc2  =  exp(-xinflam*TIME-tinflamadmin2) 

lconc3  =  exp(-xinflam*TIME-tinflamadmin3) 

lconc4  =  exp(-xinflam*TIME-tinflamadmin4) 

lconc5  =  exp(-xinflam*TIME-tinflamadmin5) 

leone  =  Iconcl  +lconc2+lconc3+lconc4+lconc5 

leff  =1.1 

xinflam=.1 

tinflamadmin1=0 

tinflamadmin2=24 

tinflamadmin3=48 

tinflamadmin4=72 

tinflamadmin5=96 

ldose=1 

Ca  =  1 


GSH  =  1 .25*exp(-.007*Systemic) 

k_split  =  0.004 
k_bind  =  0.04 
k_repair  =  0.003 
k  kill  =  0.003 
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Section  14. 


Ocular  Source  Code 


METHOD  RK4 

STARTTIME  =  0 
STOPTIME=1400 
DT  =  0.02 

{Exposure  calculations} 

Vapor  =100  ;Average  concentration  in  mg/m3 

Liquid  =  0  ;Average  concentration  in  mg/m2 

Tmask  =  1/60 

Tmopp=  1/60 

CF  =  1.2 

Sexposed=.93 

BreathingRate=0.9 

d/dt(Ocular)  =  IF  TIME  <  Tmask  THEN  Vapor  -  k_bind*Ocular*NN  else  -  k_bind*Ocular*NN 
init  Ocular  =  0 

;Equations  for  Eye  Injury 

d/dt(NN)  =  k_split*(1-NN)*NN  -  k_bind*Ocular*NN+k_repair*PARP*NAD*NI 
init  NN=1 .0 

Eyelnjury  =  1-NN 

d/dt(NI)  =  k_bind*Ocular*NN  -  k_repair*PARP*NAD*NI  -  k_kill*Protease*lnflam*Ca*NI  -  k_split*(1 -N l)*N I 
init  Nl  =  0.0 

d/dt(NKN)  =  k_kill*Protease*lnflam*Ca*NI 
init  NKN  =  0.0 

d/dt(NKA)  =  k_split*(1-NI)*NI 
init  NKA  =  0.0 

Totallnjury  =  NKA+NKN 

A  =1 .5*NAD 
PARP  =5.9088*NI 

NAD  =1 .0827*EXP(-.51 1*PARP) 

Protease  =  35.6*exp(-3.4*A) 

Inflam  =  ,9811*exp(0.0638*NI) 


Ca  =  1 


ksplit  =  0.004 
k_bind  =  0.4 
k_repair  =  10 
k  kill  =  5.8 
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Section  15. 


Inhalation  Source  Code 


METHOD  RK4 

STARTTIME  =  0 
STOPTIME=1400 
DT  =  0.02 

{Exposure  calculations} 

Vapor  =  100  ;Average  concentration  in  mg/m3 

Liquid  =  0  ;Average  concentration  in  mg/m2 

Tmask  =  10/60 

Tmopp=  2/60 

CF  =  1.2 

Sexposed=.93 

BreathingRate=0.9 

d/dt(lnhaled)  =  IF  TIME  <  Tmask  THEN  Vapor*BreathingRate  -  NACreduction  *k_bind*lnhaled*NN  else  - 
NACreduction  *k_bind*lnhaled*NN 
init  Inhaled  =  0 

NACreduction  =  1 

;Equations  for  Inhalation  Injury 

d/dt(NN)  =  k_split*(1-NN)*NN  -  k_bind*lnhaled*NN+k_repair*PARP*NAD*NI 
init  NN=1.0 

Lunglnjury  =  1-NN 

d/dt(NI)  =  k_bind*lnhaled*NN  -  k_repair*PARP*NAD*NI  -  k_kill*Protease*lnflam*Ca‘NI  -  k_split*(1-NI)*NI 
init  Nl  =  0.0 

d/dt(NKN)  =  k_kiirProtease*lnflam*Ca*NI 
init  NKN  =  0.0 

d/dt(NKA)  =  k_split*(1-NI)*NI 
init  NKA  =  0.0 

Totallnjury  =  NKA+NKN 

A  =1 .5*NAD 
PARP  =5.9088*NI 

NAD  =1.0827*EXP(-.511*PARP) 

Protease  =  35.6*exp(-3.4*A) 

Inflam  =  .981 1*exp(0.0638*NI) 


Ca  =  1 


ksplit  =  0.004 
k  bind  =  0.35 
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krepair  =  10 
k  kill  =  2 
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Appendix  A. 


Inhalation  and  Ocular  Model 


This  Appendix  preserves  the  work  done  for  inhalation  and  ocular  injury  modeling. 

A.l  Exposures 
A.  1 . 1  Inhalation 

Inhalation  dosage  is  defined  as  the  integral  of  the  vapor  concentration  (Q  from  the  start  of 
exposure  (to),  until  the  time  a  mask  (tmask)  is  donned  or  until  the  end  of  the  exposure  (t£), 
whichever  is  lower.  The  inhalation  dose  is  the  dosage,  times  the  breathing  rate  (b).  Therefore,  the 
inhalation  dose  is  defined  as: 

( j  f^mask  ( A-l ) 

cdt 

Figure  A-l  shows  the  inhaled  dose  as  a  function  of  time,  for  an  individual  (breathing  at  a  rate 
of  0.015  m  /min)  exposed  to  1  mg/m  of  vesicant  vapor,  until  the  mask  is  donned  at  2  minutes. 

As  you  can  see,  exposure  to  a  constant  vapor  results  in  a  constant  increase  in  inhaled  dose,  until 
the  mask  is  donned. 


0.035 

Inhaled  Dose  Over  Time 

0.025 

73 

£  0.02 

O) 

_ / _ 
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/ _ 

/ 

_ L _ 

/ 

c 

2  4  6  8  10 

Time  (Minutes) 

Figure  A-l.  Inhaled  Dose  Over  Time  for  a  Constant  Exposure  and  Breathing  Rate 
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Note  that  inhaled  liquid  droplets  are  not  included  in  the  model  at  this  time.  An  inhaled  liquid 
droplet  model  for  HD  has  been  developed  for  the  Joint  Effects  Model  (JEM)  and  could  be  used 
in  the  front-end  for  the  final  implementation  of  the  GYM. 


A.  1 .2  Ocular 

Ocular  dose  is  simply  the  integral  of  the  vapor  concentration  (Q  from  the  start  of  exposure 
(t0),  until  the  time  a  mask  (tmask)  is  donned  or  until  the  end  of  the  exposure  (tf),  whichever  is 
lower.  There  is  no  modification  for  breathing  rate  when  considering  ocular  exposures.  Therefore, 
the  ocular  dose  is  defined  by; 


^  r  Tmask  (A-2) 

TtD°‘l  cdt 


Figure  A-2  shows  the  ocular  dose  as  a  function  of  time,  for  an  individual  exposed  to  1  mg/m 
of  vesicant  vapor,  until  the  mask  is  donned  at  2  minutes.  As  you  can  see,  exposure  to  a  constant 
vapor  results  in  a  constant  increase  in  inhaled  dose,  until  the  mask  is  donned. 


Ocular  Dose  Over  Time 

2.5  | — 
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ID  1.5 
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Time  (Minutes) 


Figure  A-2.  Ocular  Dose  Over  Time  for  a  Constant  Exposure  and  Breathing  Rate 

At  the  time  of  this  report,  we  have  no  model  for  liquid  ocular  exposure.  For  the  foreseeable 
future,  ocular  doses  will  only  be  from  vapor  exposures. 
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A.2  Using  Ocular  and  Inhaled  Doses 


Inhaled  dose  is  used  by  calculating  the  dose  as  a  function  of  time  and  subtracting  the  rate  at 
which  the  vesicant  attacks  healthy  cells.  This  modifies  equation  (  A-3  )as  shown  below: 

^  Amask  (A-3) 

—  Di=j  C  dt  *  b  —  kbind  *  NNji  *  Dj 


Note  some  changes  in  variable  names  in  equation  (  A-3  )when  compared  to  the  original  HD 
document.  For  clarification,  we  specify  the  proportion  of  normal,  healthy  cells  in  the  lung  as 
Nnj.  The  other  tenn  introduced  above,  kbind,  is  the  rate  at  which  the  vesicant  binds  with,  and 
fragments,  the  DNA  in  nonnal,  healthy  cells.  We  fit  this  parameter  in  the  original  HD  model,  but 
we  will  be  revisiting  all  fits  in  the  GVM. 

Ocular  dose  is  used  by  calculating  the  dose  as  a  function  of  time  and  subtracting  the  rate  at 
which  the  vesicant  attacks  healthy  cells.  This  modifies  equation  (  A-3  )  as  shown  below: 


1 


Tmask 


C  dt 


kbind  *  ^N,o  *  Do 


(A-4) 


Once  again,  note  some  changes  in  variable  names  in  equation  (  A-4  ),  when  compared  to  the 
original  HD  document.  For  clarification,  we  specify  the  proportion  of  nonnal,  healthy  cells  in  the 
eye  as  NNj0. 


A.3  Injury  Model 

For  ocular  and  inhalation  exposures  we  are  using  the  same  injury  mechanism  used  for  the 
HD  model:  DNA  alkalization  and  its  effects.  We  did  make  two  significant  changes,  however,  as 
follows: 

1 .  We  included  the  depletion  of  ATP  and  correlated  that  value  with  a  pathway  to  necrotic  cell 
death. 

2.  It  is  now  the  depletion  of  ATP  in  the  cell  that  signals  protease,  kills  the  cells,  and  begins  the 
vesication  process.  In  the  HD  model,  it  was  the  depletion  of  NAD+  that  directly  signaled 
protease  release. 

These  changes  were  made  to  better  fit  the  revised  injury  mechanism  of  vesicants,  as 
presented  in  the  previous  section  and  to  better  incorporate  the  beneficial  effects  of 
countermeasures.  Looking  at  the  figure  below,  it  is  apparent  that  the  new  framework  would 
better  account  for  PARP  inhibitors.  PARP  inhibitor  would  slow  the  depletion  of  NAD+  and 
ATP,  thereby  slowing  the  release  of  protease.  This  would  result  in  the  main  path  to  necrosis 
occurring  during  the  attempted  mitosis  of  injured  cells.  In  the  percutaneous  model,  which  we 
present  below,  attempted  mitosis  of  injured  cells  leads  to  apoptotic  cell  death,  which  is  more  in 
agreement  with  the  known  injury  mechanism  presented  in  the  previous  section.  We  believe  it  is 
beneficial  to  achieve  that  level  of  fidelity  for  the  percutaneous  exposures,  due  to  the  wider  dose 
range  of  toxic  skin  exposures  and  the  fact  that  blistering,  which  results  from  necrotic  death,  is  a 
significant  endpoint  for  skin  exposures.  Additionally,  we  have  data  that  supports  a 
necrosis/apoptosis  model  for  skin  injuries,  but  not  for  ocular  and  lung  injuries. 
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Eye  or  Lung 


Combined 

Exposure 


Intracellular 


Figure  A-3.  GYM  Model  Framework  for  Ocular  and  Inhalation  Exposures 


A.4  Proportion  of  Healthy  Cells 

Given  the  above  statements,  we  can  draw  the  system  diagram  for  Nn,l  and  Nn,o,  shown 
inFigure  A-4. 


Figure  A-4.  Proportion  of  Healthy  Cells  in  the  Lungs  (Nn,l)  and  Eyes  (Nn,q) 


The  equations  that  describe  the  above  diagram  are: 

d  (A-5 ) 

=  ksput  *  (l  —  Nn  l )  *  Nn  l  —  kbind  *  Di  *  Nn  l  +  krepair  *  PARPL  *  NADL  *  NI  L 
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(A-6) 


foNN,o  ~  k spin  *  (l  Nn  0)  *  Nn  0  kbind  *  D0  *  Nn  o  +  krepair  *  PARP0  *  NAD0  *  NI0 
Where: 

1 .  D;  and  D0  are  the  inhaled  and  ocular  doses,  respectively,  given  in  milligrams  and  described  in 
Section  5, 

2.  kgpiit  is  the  rate  of  mitotic  division,  given  in  min' , 

3.  kbmd  is  the  rate  vesicant  binds  with  the  DNA  in  healthy,  normal  cells  in  min' , 

4.  krepair  is  the  rate  of  DNA  repair,  given  in  min' ,  and 

5.  NADl  and  NAD0  are  the  relative  concentrations  of  NAD+  and  PARP  in  the  lungs  and  eyes, 
respectively,  and  are  defined  by  the  equations  in  Section  6.2.1. 

A.5  Proportion  of  Injured  Cells 

Given  the  above  statements,  we  can  draw  the  system  diagram  for  Ni,l  and  Ni,o,  shown  in 
Figure  A-5. 


Figure  A-5.  Proportion  of  Injured  Cells  in  the  Lungs  (NN;L)  and  Eyes  (NNi0) 


The  equations  that  describe  the  above  diagram  are: 

d  (A-7) 

—  kbind  *  Nnl  *  Dj  krepair  *  N/i  *  NAD  *  PARP  kkm  *  Prot  *  Ca  *  Inf  *  NIL 

~  k split  *  (l  —  N,  L)  *  Ni  l 

d  (A-8) 

—  kbind  *  Nn  0  *  D0  krepair  *  NI0  *  NAD  *  PARP  kkiu  *  Prot  *  Ca  *  Inf 

*  N. i,o  ~  k sPut  *  (l  _  Ni  o)  *  Nj  0 

Where: 
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1 .  D;  and  D0  are  the  inhaled  and  ocular  doses,  respectively,  given  in  milligrams, 

2.  kgpiit  is  the  rate  of  mitotic  division,  given  in  min'1, 

3.  kbmd  is  the  rate  vesicant  binds  with  the  DNA  in  healthy,  normal  cells  in  min'1, 

4.  krepair  is  the  rate  of  DNA  repair,  given  in  min'1, 

5.  kjdu  is  the  rate  of  DNA  death,  given  in  min'1,  and 

6.  NAD,  PARP,  Prot,  Ca,  and  Inf  are  the  relative  concentrations  of  NAD+,  PARP,  protease, 
calcium  ions,  and  inflammatory  cells,  respectively,  as  defined  in  section  6.2. 

A.6  Proportion  of  Necrotic  Cells 

We  can  draw  a  system  diagram  for  necrotic  cell  death  in  the  lungs  and  eyes,  as  shown  in 
Figure  A-6. 


Figure  A-6.  Proportion  of  Necrotic  Cells  in  the  Lungs  (Nk,l)  and  Eyes  (NK,o) 


The  equations  that  describes  the  above  diagram  are: 


d  .  . 

—  iVKL  =  ksplit  *  (1  -  NjL J  *  ni,l  +  kkui  *  Prot  *Ca*  Inflam  *  NIL 

d  ,  . 

—  =  ksplit  *  (1  -  NIi0)  *  Nj  0  +  kkill  *  Prot  *Ca*  Inf  lam  *  N,i0 
Where: 

1.  kSpiit  is  the  rate  of  mitotic  division,  given  in  min'1,  and 

2.  kkm  is  the  rate  of  necrotic  cell  death,  given  in  min'1 


(A-9) 
(  A-10  ) 


A.7  The  Repair  Rate,  Kill  Rate,  and  HD  Binding  Rate  for  Ocular  Exposures 

For  ocular  exposures,  we  assumed  that  the  mitotic  division  rate  is  still  0.004  hr'1,  since  we 
found  no  data  that  supported  changing  the  rate.  We  assumed  that  the  other  values  could  change, 
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due  to  differences  in  the  cellular  structure  of  the  eye.  Holding  that  value  constant,  we  fit  a 
10mg/mJ  exposure  for  ten  minutes  to  the  highest  ocular  dose  band  in  AMedP-8  (C).  At  that  dose 
band,  ocular  injury  peaks  at  1 1  hours  and  recovers  at  672  hours.  The  curves  corresponding  to  our 
lit  can  be  seen  in  Figure  A-7.  The  values  that  correspond  to  the  curve  are  tabulated  in  Table  A-l. 


Figure  A-7.  Ocular  Model  Output  for  HD  Compared  to  AMedP-8  (C) 


Table  A-l.  HD  Ocular  Values  for  kbind,  krepair,  and  kkni 


kbind  (hr  ) 

krepair  (hr 
!) 

kkni  (hr' 
l) 

0.4 

10.0 

5.8 

Note  the  higher  values  for  kbind,  krepair,  and  ktai,  compared  to  the  percutaneous  exposure.  This 
is  due  to  the  earlier  peak  time  (1 1  hours  vs.  24  hours)  and  the  lower  dosage  (100  mg  vs.  180  mg) 
used  in  the  fit. 

A.8  The  Repair  Rate,  Kill  Rate,  and  HD  Binding  Rate  for  Inhalation  Exposures 

For  inhalation  exposures  we  also  assumed  that  the  mitotic  division  rate  is  still  0.004  hr'1, 
since  we  found  no  data  that  supported  changing  the  rate.  We  assumed  that  the  other  values  could 
change,  due  to  differences  in  the  cellular  structure  of  the  lung.  Holding  that  value  constant,  we  fit 
a  25  mg/m3  exposure  for  10  minutes  to  the  150-250  mg/m3  inhalation  dose  band  in  AMedP-8 
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(C).  At  that  dose  band,  lung  injury  peaks  at  36  hours  and  recovers  at  1008  hours.  The  curves 
corresponding  to  our  fit  can  be  seen  in  Figure  A-8.  The  values  that  correspond  to  the  curve  are 
tabulated  in  Table  A-2. 


Table  A-2  -  HD  Inhalation  Values  for  kbind,  krepair,  and  kkiu 


kbind  (hr  ) 

krepair  (hr  ) 

kkiii  (hr1) 

0.35 

10.0 

2 

Note  that  the  values  for  khind  and  krepair  are  equal  to  the  ocular  values  and  kkiii  is  the  only 
different  value.  This  is  due  to  the  slower  recover  time  (1008  hours  vs.  672  hours)  used  in  the  fit. 
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